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(54) Interpolating operation method and apparatus for image signals 

from those of the original image signal components Yij, 



(57) An original image signal, which represents an 
original image and is made up of a series of original 
image signal components Yij, is obtained. Correspond- 
ing interpolation coefficients Bij and Cij, which are set 
for each of the original image signal components Yij, in 
two different interpolating functions f and g, represented 
by Formulas (1) and (2), for obtaining two interpolation 
images having different levels of sharpness are linearly 
combined with each other by using Formula (3). A new 
interpolation coefficient Aij is thereby obtained. An inter- 
polating operation is carried out on the original image 
signal components Yij with an interpolating function h ( 
which has the new interpolation coefficient Aij and is 
represented by Formula (4). An interpolation image sig- 
nal is thereby obtained, which is made up. of a series of 
image signal components occurring at intervals different 

F I 



f=LBij • Yij 
g=LCij • Yij 
Aij=(1-a)Bij+aCij 
h=£Aij * Yij 



(1) 
(2) 
(3) 
(4) 



in which i=1, 2, ... ( and j=1, 2, 
The coefficient a in Formula (3) is set to be one of real 
numbers including a range smaller than 0 and/or a 
range larger than 1 . 
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Description 

BACKGROUND OF THE INVENTION 

5 Field of the Invention 

This invention relates to an interpolating operation method and apparatus for an image signal. 
Description of the Prior Art 

10 

Techniques for photoelectrical^ reading out an image, which has been recorded on a photographic film, in order to 
obtain an image signal, carrying out appropriate image processing on the image signal, and then reproducing a visible 
image by use of the processed image signal have heretofore been known in various fields. 

Also, it has been proposed to use stimulable phosphors in radiation image recording and reproducing systems, 

is Specifically, a radiation image of an object, such as a human body, is recorded on a sheet provided with a layer of the 
stimulable phosphor (hereinafter referred to as a stimulable phosphor sheet). The stimulable phosphor sheet, on which 
the radiation image has been stored, is then exposed to stimulating rays, such as a laser beam, which cause it to emit 
light in proportion to the amount of energy stored thereon during its exposure to the radiation. The light emitted by the 
stimulable phosphor sheet, upon stimulation thereof, is photoelectrical^ detected and converted into an electric image 

20 signal. The image signal is then processed and used for the reproduction of the radiation image of the object as a visible 
image on a recording material, such as photographic material, or on a display device, such as a cathode ray tube (CRT) 
display device. Radiation image recording and reproducing systems, which use stimulable phosphor sheets, are advan- 
tageous over conventional radiography using silver halide photographic materials, in that images can be recorded even 
when the energy intensity of the radiation, to which the stimulable phosphor sheet is exposed, varies over a wide range. 

25 In image recording and reproducing systems, in which an image signal is obtained in the manner described above 
and a visible is reproduced from the image signal, in cases where the region of interest in the visible image is to be 
viewed in more detail, the region of interest is often enlarged and reproduced. In such cases, if the enlargement of the 
image size is carried out such that the number of the image signal components of the image signal representing the 
enlarged image may be identical with the number of the image signal components of the original image signal repre- 

30 senting the original image, the sharpness of the enlarged image will be recognized to be lower than the sharpness of 
the original image due to the visual characteristics of persons. Therefore, if the image is merely enlarged and repro- 
duced, an enlarged image having a high sharpness cannot be obtained, and the details of the image cannot be viewed 
accurately. 

In order for the aforesaid problems to be eliminated, a predetermined interpolating operation may be carried out on 
35 the original image signal, which has been obtained by reading out an original image, and an interpolation image signal, 
which is a secondary image signal and is made up of a number of image signal components different from that of the 
original image signal, may thereby be formed. Specifically, in cases where an enlarged image is to be reproduced, an 
interpolation image signal, which is made up of a number of image signal components larger than that of the original 
image signal, may be formed from the interpolating operation. A visible image may then be reproduced from the inter- 
ne polation image signal. In this manner, the sharpness of the enlarged image can be prevented from becoming low. 

As the interpolating operation methods for carrying out interpolating operations on image signals, various methods 
have heretofore been proposed. Among such methods, the method using third-order spline interpolating functions is 
popular. With the interpolating operation method using the third-order spline interpolating functions, digital original 
image signal components (Y^ in each section are connected by a third-order function {f k }, and the value of f k corre- 
45 sponding to a position, at which an interpolation point is set, (i.e., a setting position in each section) is taken as the value 
of the interpolated image signal component 

The interpolating operations, which pass through the original image signal in the manner described above, can 
yield an image having a comparatively high sharpness. As such interpolating operations, cubic spline interpolating 
operations, and the like, are known. How the cubic spline interpolating operations are carried out will be described here- 
so inbelow. 

Figure 4 is an explanatory graph showing how interpolated image signal components are obtained with a conven- 
tional cubic spline interpolating operation from original image signal components, which are sampled with a period of 
an equal interval and represent sampling points (picture elements) arrayed in one direction. As illustrated in Figure 4, 
the image signal components (the original image signal components), which have been detected as digital signal com- 
55 ponents from an original image and represent a series of picture elements X^, , X^ X|< +1 , X^, are respectively 
represented by Y k _ 2 , Y k _ 1f Y k , Y k+1 , Y k+2 , ... A third-order spline interpolating function is set for each of sections X^- 
^k-i • *k-i ~ x k. Xk-Xk+i • anc * *k+i -^+2- The spline interpolating functions corresponding to the respective sections are 
represented by f k . 2t fk-i. fk. f k+i. and f k+2 . The interpolating functions are the third-order functions, in which the position 
in each section serves as a variable. 
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Firstly, how the interpolating operation is carried out when a point taken for interpolation (hereinbelow referred to 
as the interpolation point) Xp falls within the section X^ X,< +1 will be described hereinbelow. The spline interpolating 
function f k corresponding to the section X k ~X k+1 is represented by Formula (29). 

5 f k (x)=A k x 3 +B k x 2 +C k x+D k (29) 

In the cubic spline interpolating operation, it is necessary that the spline interpolating function f k passes through the 
original sampling points (picture elements), and that the first-order differential coefficient of the spline interpolating func- 
tion f k is continuous between adjacent sections. Therefore, it is necessary for Formulas (7), (8), (9), and (10) to be sat- 
w isfied. 

MXk>=Yk (7) 
M* K+1 )=Y k+1 (8) 

W k+1 )=f k+1 xx k+1 ) do) 

20 In these formulas, f k ' represents the first-order differentiation (3A k x 2 +2B k x+C k ) of the function f k . 

In the strict sense, the cubic spline interpolating operation contains the continuity conditions of the second-order 
differential coefficient. However, with continuity conditions of the second-order differential coefficient, the operation for- 
mulas become complicated. Therefore, the cubic spline interpolating operation is popularly carried out in the form sim- 
plified in the manner described above. 

25 Also, in the cubic spline interpolating operation, it is necessary for the first-order differential coefficient at the picture 
element Xj< to satisfy the condition with respect to the picture elements X k . 1 and X k+1 , which are located before and after 
the picture element X k , in that the first-order differential coefficient at the picture element X|< should coincide with the 
gradient (Y k+1 i Y k . 1 )/(X k+1 -X k .-,) of the image signal components Y k . 1 and Y k+1 representing the picture elements X^ 
! and X^. Therefore, it is necessary for Formula (19) to be satisfied. 
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fkXX^Y^-Y^^X^-X^) (19) 



Also, it is necessary for the first-order differential coefficient at the picture element X^ to satisfy the condition with 
respect to the picture elements X k and X^, which are located before and after the picture element X k+1 , in that the f irst- 
35 order differential coefficient at the picture element X k+1 should coincide with the gradient (Y k+2 -Y k )/(X k+2 -X k ) of the 
image signal components Y k and Y k+2 representing the picture elements Xj< and X k+2 - Therefore, it is necessary for For- 
mula (20) to be satisfied. 



f k XX k+1 )=(Y k+2 -Y k )/(X k+2 -X k ) (20) 

It is herein assumed that the interval (i.e., the lattice interval) of each of sections X k . 2 -X k . 1t X^-X^ X k - X k+1 , and 
X k+ i-X k+2 is equal to 1 . and the position of the interpolation point Xp, which is taken from the picture element X* toward 
the picture element Xk +1 , is represented by t (Osts 1). In such cases, from Formulas (29), (7), (8). (9), (10), (19), and 
(29), the formulas shown below obtain. 



f k (0)=D k =Y k 
f k (1)=A k +B k+ C k+ D k =Y k+1 

50 WO^Ck-O^ 

f k (1)=3A k+ 2B k+ C k =(Y k+2 -Y k )/2 

Therefore, the formulas shown below obtain. 

AK^K+a-SY^+aY.-YK.!)* 

B k -(-Y lHa+ 4Y lw1 -5Y k rfY k . 1 >/2 
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As described above, the variable conversion of X=t is carried out, and therefore the spline interpolating function 
f k (x) is represented by the formula shown below. 

*k<xM k (t) 

Therefore, an interpolated image signal component Y p corresponding to the interpolation point X p may be represented 
by Formula (30). 

Y P= f k(t)=A k t 3 +B k t 2 + C k t + D k (30) 
Substituting the coefficients \, B k , C k , and D k into Formula (30) yields 

Y p =((Y k+2 -3 Y k+1 + 3Y k -Y ko )/2}t 3 +{(-Y k+2 +4Y k+1 -5 Y k +2Y k _, )/2}t 2 + {(Y k+ r Y h , )/2}t + Y k 

Arranging this formula with respect to the image signal components Y k . 1( Y k , Y k+1 , and Y k+2 yields Formula (31). 

Y p={(-t 3 +2t 2 -t)/2}Y k +{(3t 3 -5t 2 +2)/2}Y k +{(-3t 3 +4t 2 +t)/2}Y k+1 +{(t 3 -t 2 )/2}Y k+2 (31) 

The coefficients for the image signal components Y k _., , Y k , Y k+1 , and Y k+2 are referred to as the interpolation coef- 
ficients cvl c k , C(< + 1 , and q^. Specifically, the interpolation coefficients c^, c^ c k+1 , and c^, which respectively cor- 
respond to the image signal components Y k . 1( Y k , Y k+1 , and Y k+2 in Formula (31), may be represented by the formulas 
shown below. 

c^t-t 3 -^ 2 -^ 
c k =(3t 3 -5t 2 +2)/2 
c k+1 =(-3t 3 +4t 2 +t)/2 
c k+2 =(t 3 -t 2 )/2 

The operations described above are repeated for the sections X^-X^, X^-X^ X k -X k+1( and X k+1 ~X k+2 . In this 
manner, an interpolation image signal can be obtained, which is made up of image signal components occurring at 
intervals different from those of the image signal components of the entire original image signal. 

As described above, in the cubic spline interpolating operation, it is necessary that the spline interpolating function 
passes through the original sampling points (picture elements), and that the first-order differential coefficient of the 
spline interpolating function is continuous between adjacent sections. With the interpolating function for the cubic spline 
interpolating operation, the interpolation image signal for use in the reproduction of a secondary image (i.e., the image 
obtained from the interpolating operation), which has a comparatively high sharpness, is obtained. On the other hand, 
as for a portion in the original image, at which the change in density is gentle, the interpolating operation should pref- 
erably be carried out such that a secondary image, in which the sharpness is comparatively low and which is smooth, 
may be obtained. As the interpolating function for obtaining the interpolation image signal representing the secondary 
image, in which the sharpness is comparatively low and which is smooth, for example, a B spline interpolating operation 
function is known. In the B spline interpolating operation, the spline interpolating function need not pass through the 
original sampling points (picture elements), and it is necessary that the first-order differential coefficient and the sec- 
ond-order differential coefficient {represented by f '(X)} of the spline interpolating function are continuous between adja- 
cent sections. 

Specifically, in Formula (29), 

f k (x)=A k x 3 +B k x 2 +C k x+D k (29) 
the conditions shown below should be satisfied. 



V(X k )=W(X k ) 
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V(x k+1 ')=W(X k+ i) < 10 > 
VCXkMk-iW (32) 

5 fK"(X k+1 )=f k+n "(X k+1 ) (33) 

Also, it is necessary for the first-order differential coefficient at the picture element X k to satisfy the condition with 
respect to the picture elements X^ and X k+1 , which are located before and after the picture element X k , in that the first- 
order differential coefficient at the picture element X k should coincide with the gradient (Y k+1 -Y ^/(X k+1 -X M ) of the 
10 image signal components Y k _-| and Y k+1 representing the picture elements X k ^ and X k+1 . Therefore, it is necessary for 
Formula (19) to be satisfied. 

fk'(X k )=(Y k+l -Y k . l )/(X k+1 -X k . l ) (19) 



Further, it is necessary for the first-order differential coefficient at the picture element X k+1 to satisfy the condition 
with respect to the picture elements X k and X k+2 , which are located before and after the picture element X k+1 , in that 
the first-order differential coefficient at the picture element X,< +1 should coincide with the gradient (Y k+2 -Y k )/(X k+2 -X k ) 
of the image signal components Y k and Y k+2 representing the picture elements X* and X^. Therefore, it is necessary 
for Formula (20) to be satisfied. 

fK(X k+n )=(Y k+2 -Y k )/(X k+2 -X k ) (20) 



In general, the function f(X) may be approximately represented by Formula (34). 
25 f(X)=f(0)+f(0)X+{f\0)/2}X z (34) 

It is herein assumed that the interval (i.e.. the lattice interval) of each of sections X^-X^.,, X^-X^ Xk- X^,, and 
Xk +1 -Xk+2 is equal to 1 , and the position of the interpolation point Xp, which is taken from the picture element Xk toward 
the picture element X k+1 . is represented by t (Osui 1). In such cases, from Formulas (29), (9), (10), (19), (20), (30), (32), 
30 (33), and (34). the formulas shown below obtain. 

f k '(1 )=3A k +2B k+ C k =(Y k+2 -Y k )/2 

35 

f K "(0)=Y k+1 -2Y k+ Y k . 1 =2B 
Therefore, the formulas shown below obtain. 
« A^O^-SY^+SYk-Y^/e 

B k =(Y k+1 -2Y k+ Y k .,)/2 

45 

Since D k is unknown, it is represented by the formula * 

D k -(D , Y k+2+ D 2 Y k+1 +D 3 Y k+ D 4 Y M )/6 

so As described above, the variable conversion of X=t is carried out, and therefore the spline interpolating function 
f k (x) is represented by the formula shown below. 

f k (x)=f k (t) 

55 Therefore, 

f k (t)=((Y k+2 -3Y k+1 +3Y k -Y M )/6}t 3 + {(Y k+1 -2Y k+ Y M )/2}t 2 
+{(Y k+1 -Y m )/2}t +(D 1 Y k+2+ D 2 Y k+1 +D 3 Y k+ D 4 Y M )/6 
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Arranging this formula with respect to the image signal components Y k . 1f Y k . Y k+1 , and Y k+2 yields Formula (35). 

^W^C-t^at^St+D^/ejY,., + {(3t 3 -6t 2 + D 3 )/6}Y k +{(-3t 3 + 3t 2 + 3t + D 2 )/6}Y k+1 + {(t 3 + D .J/ejY^ (35) 
If t is set to be t=1 , the formula shown below will obtain. 

f k (1 H(D 4 -1 )/6}Y M + {(D 3 -3)/6}Y k +{(D 2+ 3)/6}Y k+1 + {(D , +1 )/6}Y k+2 
As for the section X^-X^, Formula (35) may be rewritten as Formula (36) 

f k+ i(t)={(-t 3 + 3t 2 -3t + D 4 )/6}Y k +{(3t 3 -6t 2 + D 3 )/6}Y k+1 + {(-3t 3 + 3t 2 + 3t + D 2 )/6}Y k+2 + {(t 3 + D , )/6}Y k+3 (36) 
If t is set to be t=0, the formula shown below will obtain. 

f k+1 (0)=(D 4 /6)Y k+ (D 3 /6)Y k+1 + (D 2 /6)Y k+2+ (D ^Y^ 

From the continuity condition {f k (l )=f k+1 (0 ) } and the condition in that the coefficients corresponding to the respec- 
tive original image signal components are equal to each other, D 4 -1=0, D 3 -3=D 4 , D 2 +3=D 3 , D 1+ 1-D 2 , and D 1= 0. 



D k=(Yk + i+4Y k+ Y k . 1 )/6 



Accordingly, Formula (37) obtains. 



Y p «f k (t)={(-t 3 +3t 2 -3t+1)/6}Y k . 1 +{(3t 3 -6t 2 +4)/6}Y k -K(-3t 3 + 3t 2 + 3t+1)/6}Y k+1 +(t 3 /6)Y k+2 (37) 

Therefore, the interpolation coefficients b^ , k^ +1 , and b^, which respectively correspond to the image signal 
components Y M , Y k , Y k+1 , and Y k+2 , may be represented by the formulas shown below. 



30 K _/*3.o*2 



b k-i=(-t +3t -3t+1)/6 
b k =(3t 3 -6t 2 +4)/6 
b k+1 =(-3t 3 +3t 2 +3t+1)/6 



b k+2 =t 3 /6 
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man™ Tnll^ ^ V tor th6 SeCt, '° ns **-*~X*+ *k~X k+1 . and Xk +1 ~Xk +2 . In this 

manner, an interpolation image signal can be obtained, which is made up of image signal components occurring at 
intervals different from those of the image signal components of the entire original image signal 

In this manner, in cases where a secondary image (an interpolation image) having a high sharpness is to be repro- 
duced for example, the cubic spline interpolating operation may be used. In cases where a secondary image which 
has a low sharpness and is smooth, is to be reproduced, for example, the B spline interpolating operation may be used 

in U.5.R 5.048. 105. the applicant proposed an interpolating operation method for an image signal, with which the 
sharpness of an interpolation image can be adjusted finely by. for example, weighting the corresponding coefficients of 
two mterpolating functions, that provide different levels of sharpness, in accordance with a desired sharpness of the 
interpolate image, and adding the weighted coefficients to each other. With the proposed method, for example in 
cases where the cubic spline interpolating operation function and the B spline interpolating operation function are 
employed as the two interpolating functions, that provide different levels of sharpness, the interpolation coefficients c* 
I! °£ ^S 1 ' r f ° r CUb,C Sp ' ine inter P° latin 9 operation and the interpolation coefficients bu , , bk. b^, . and b^ 9 
for the B spline interpolating operation, which coefficients correspond to each other and are set for thrrel£ective orig- 
inal .mage signal components Y^. Y k , Y k+1 , and Y k+2 . are weighted and added to each other. By alteration of the 
we.ghtng ratio (coefficient) a. it is possible to obtain a secondary image having a desired level of sharpness, which is 
intermediate in the range from the highest sharpness to the lowest smooth sharpness 

Specifically, in cases where the interpolation coefficients for the cubic spline interpolating operation are repre- 
sented by q,., . ck. c^,. and c*. z , and the interpolation coefficients for the B spline interpolating operation are repre- 
sented by bk., , bk. bk +1 . and bk +2 . weighted interpolation coefficients a k .„ a*, a**,, and a^ are set as shown below. 

a k-1=( 1 -«)Ck-l +ctb k-1 
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a k =(1.a)c k +alD k 



a k+1 =(1-a)c k+1 +ctb 



k+1 



a k + 2=( 1 - a ) C k + 2 +ab 



k+2 



wherein 0^a^1 . 

In accordance with the thus obtained new interpolation coefficients a k . 1( a k , a k+1 , and a k+2 > an interpolated image 
signal component Y p is calculated with Formula (38). 



An actual image is composed of the picture elements arrayed in two-dimensional directions. Therefore, the interpo- 
lation coefficient a k is represented as the interpolation coefficient Bij or Cij corresponding to each of two different direc- 
tions (an i direction and a j direction) of an array of picture elements in the image. 

With respect to the sharpness of the interpolation image, a wider variety of sharpness levels are often desired. For 
example, it is often desired to obtain an interpolation image having a sharpness higher than the sharpness of the inter- 
polation image, which is obtained from the cubic spline interpolating operation alone. Also, it is often desired to obtain 
an interpolation image having a sharpness smoother than the sharpness of the interpolation image, which is obtained 
from the B spline interpolating operation alone. 

However, with the interpolating operation method for an image signal, which is disclosed in U.S. P. 5,048,105, in 
cases where the cubic spline interpolating operation and the B spline interpolating operation are employed, the adjust- 
ment of sharpness can be carried out only within the range of the sharpness, which corresponds to the sharpest image 
obtained with the cubic spline interpolating operation, to the sharpness, which corresponds to the smoothest image 
obtained with the B spline interpolating operation. Therefore, the disclosed method cannot satisfy the demand for a 
wide variety of sharpness levels. 

SUMMARY OF THE INVENTION 

The primary object of the present invention is to provide an interpolating operation method for an image signal, 
wherein the sharpness of an interpolation image obtained by linear combination of two kinds of interpolating operations, 
that provide different levels of sharpness, is capable of being adjusted with a high degree of freedom. 

Another object of the present invention is to provide an apparatus for carrying out the interpolating operation 
method for an image signal. 

A further object of the present invention is to provide a spline interpolating operation method for an image signal, 
wherein the degree of freedom of sharpness adjustment is kept high only with a single kind of spline interpolating oper- 
ation. 

A still further object of the present invention is to provide an apparatus for carrying out the spline interpolating oper- 
ation method for an image signal. 

The present invention provides an interpolating operation method for an image signal, comprising the steps of: 

i) obtaining an original image signal, which represents an original image and is made up of a series of original 
image signal components Yij, 

ii) linearly combining interpolation coefficients Bij and Cij, which correspond to each other and are set for each of 
the original image signal components Yij, in two different interpolating functions f and g for obtaining two interpola- 
tion images having different levels of sharpness, which functions are represented by Formulas (1) and (2), t the lin- 
ear combination being carried out with Formula (3), a new interpolation coefficient Aij being obtained from the linear 
combination, and 

iii) carrying out an interpolating operation on the original image signal components Yij by using an interpolating 
function h having the new interpolation coefficient Aij, which function is represented by Formula (4), an interpolation 
image signal being obtained from the interpolating operation, the interpolation image signal being made up of a 
series of image signal components, which occur at intervals different from those of the original image signal com- 
ponents Yij, 



Y p =a k-1 Y k-1 +a k Y k +a k+1 Y k+1 +a k+2 Y k+2 



(38) 



f=LBij • Yij 



(1) 



g=£Cij - Yij 



(2) 



Aij=(1-a)Bij+aCij 



(3) 
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h=LAij • Yij . (4) 

in which i=1 , 2 and j=1 , 2 

wherein the coefficient a in Formula (3) is set to be one of real numbers including a range smaller than 0 
and/or a range larger than 1 . 

In the interpolating operation method for an image signal in accordance with the present invention, as for the two 
different interpolating functions for obtaining two interpolation images having different levels of sharpness, a B spline 
interpolating operation function should preferably be employed as the interpolating function corresponding to an image 
having a comparatively low sharpness, and a cubic spline interpolating operation function should preferably be 
employed as the interpolating function corresponding to an image having a sharpness higher than the comparatively 
low sharpness. This is because, in cases where the B spline interpolating operation function and the cubic spline inter- 
polating operation function are combined with each other, the first-order differential coefficient becomes continuous. 

However, the interpolating operation method for an image signal in accordance with the present invention is not lim- 
ited to the combination of the B spline interpolating operation function and the cubic spline interpolating operation func- 
tion. For example, various interpolating operation functions, such as the B spline interpolating operation function, the 
cubic spline interpolating operation function, a linear interpolating function, and a Lagrangean interpolating operation 
function, may be used, and an arbitrary combination of two of these functions may be employed. 

Each of the term "interpolation coefficient Bij** and the term "interpolation coefficient Cij" as used herein means the 
interpolation coefficient corresponding to each of two different directions (an i direction and a j direction) of an array of 
picture elements in the image. (Each of the interpolation coefficients Bij and Cij corresponds to the coefficient, by which 
each original image signal component is multiplied, e.g. each of the interpolation coefficients a*.-,, a k , a k+1l and a k+2 
described above for the conventional technique.) 

The interpolating operation method for an image signal in accordance with the present invention may be modified 
such that: 

a) a plurality of first look-up tables are prepared, each of the first look-up tables setting the frequencies and the cor- 
responding responses with respect to one of the two interpolating functions and for each of a plurality of different 
image size enlargement scale factors, 

b) a plurality of second look-up tables are prepared, each of the second look-up tables setting the frequencies arid 
the corresponding responses R 2 with respect to the other interpolating function and for each of a plurality of differ- 
ent image size enlargement scale factors, 

c) a response of the one interpolating function with respect to an image size enlargement scale factor, which is 
desired for an interpolation image, and a response R 2 of the other interpolating function with respect to the image 
size enlargement scale factor, which is desired for the interpolation image, are calculated by making reference to a 
first look-up table and a second look-up table, which correspond to the image size enlargement scale factor desired 
for the interpolation image, and 

d) the coefficient a is calculated from an operation carried out with Formula (5): 

a=(R-R 1 )/(R 2 -R 1 ) {5) 

in accordance with a response R desired for the interpolation image, the calculated response R 1 of the one inter- 
polating function, and the calculated response R 2 of the other interpolating function. 

An example of the first look-up table corresponding to a predetermined image size enlargement scale factor is illus- 
trated in Figure 5A. An example of the second look-up table corresponding to the predetermined image size enlarge- 
ment scale factor is illustrated in Figure 5B. By way of example, a frequency of 1 cycle/mm, to which the eyes of persons 
are most sensitive, is taken as a frequency of interest, and the responses Rj and R 2 at the frequency of 1 cycle/mm are 
found from the first and second look-up tables. A desired response R can then be obtained by carrying out an interpo- 
lating operation using an equation of the first order, R=aR 2 +(1-a)R 1 . Arrangement of this equation with respect to the 
coefficient a yields Formula (5). The frequency of interest is not limited to 1 cycle/mm. The other frequency may be 
selected as the frequency of interest in accordance with the kind of the image, or the like. As another alternative, two 
or more different frequencies may be selected, and the arithmetic mean value of the responses at the two or more dif- 
ferent frequencies, or the like, may be employed as the representative value. 

The extent of a change in sharpness can be easily grasped as a change in image response. Therefore, with the 
technique for directly specifying the response R as a desired value, an interpolation image having a sharpness truer to 
nature can be obtained than with the technique for specifying the value of a serving as a mere inorganic coefficient. 

In cases where the technique for directly specifying the response R is employed and the image size is to be 
enlarged, it is necessary for the image size enlargement scale factor to be specified with a certain known means. If the 



8 



EP 0 753 828 A2 



first look-up table and the second look-up table corresponding to the specified image size enlargement scale factor are 
contained in the plurality of the look-up tables having been prepared, the first look-up table and the second look-up table 
may be selected. If the first look-up table and the second look-up table corresponding to the specified image size 
enlargement scale factor are not contained in the plurality of the look-up tables having been prepared, two first look-up 

5 tables corresponding to two image size enlargement scale factors, which are closest to the specified image size 
enlargement scale factor, may be selected from the first look-up table group. A response may then be calculated by car- 
rying out an interpolating operation using an equation of the first order on the responses found from the two selected 
first look-up tables. Also, two second look-up tables corresponding to two image size enlargement scale factors, which 
are closest to the specified image size enlargement scale factor, may be selected from the second look-up table group. 

w A response may then be calculated by carrying out an interpolating operation using an equation of the first order on the 
responses found from the two selected second look-up tables. 

For example, six kinds of first look-up tables corresponding to image size enlargement scale factors of 1 .0, 1 .2 

1 .8, and 2.0 may be prepared. Also, six kinds of second look-up tables corresponding to image size enlargement scale 
factors of 1 .0, 1.2, .... 1.8, and 2.0 may be prepared. In such cases, if the specified image size enlargement scale factor 

is is, for example, 1.3, responses R^l.2) and R^l.4) may be obtained from the first look-up table for the image size 
enlargement scale factor of 1 .2 and the first look-up table for the image size enlargement scale factor of 1 A. A response 
Ri(1.3) may then be calculated by carrying out an interpolating operation using an equation of the first order, 
R^I.S^O.SR^I^J+O.SR^I^) , on the responses Ri(1.2) and R^l.4). In the same manner, the response R 2 of the 
other interpolating function may be calculated with the formula. R 2 (1 .3)=0.5R 2 (1 .2)+0.5R 2 (1 .4) . 

20 Specifically, in general, in cases where look-up tables for image size enlargement scale factors of n 1t n 2 , n 3 nj, 

and n j+1 are prepared, if the specified image size enlargement scale factor is nj (i^j^i+1 ), a response R-i(nj) may be 
obtained from the look-up table for n jt and a response Ri(n i+1 ) may be obtained from the look-up table for n M . In 
accordance with the responses Ri(nj) and Ri(n i+1 ). the response of one interpolating function may be calculated 
with the formula, R^n^s • R.|(n,)+(1-s) • R^n^) , by using the real number, s, satisfying the condition of 0^s^1. 

25 The response R 2 of the other interpolating function can be obtained in the same manner as that described above. 
The present invention also provides an interpolating operation apparatus for an image signal, wherein: 



i) an original image signal, which represents an original image and is made up of a series of original image signal 
components Yij ( is obtained, 

30 ii) interpolation coefficients Bij and Cij, which correspond to each other and are set for each of the original image 
signal components Yij, in two different interpolating functions f and g for obtaining two interpolation images having 
different levels of sharpness, which functions are represented by Formulas (1) and (2), are linearly combined with 
each other, the linear combination being carried out with Formula (3), a new interpolation coefficient Aij being 
obtained from the linear combination, and 

35 iii) an interpolating operation is carried out on the original image signal components Yij by using an interpolating 
function h having the new interpolation coefficient Aij, which function is represented by Formula (4), an interpolation 
image signal being obtained from the interpolating operation, the interpolation image signal being made up of a 
series of image signal components, which occur at intervals different from those of the original image signal com- 
ponents Yij, 

40 

f=£Bij*Yij (1) 
g=£Cij • Yij (2) 
45 Aij=(1-cc)Bij+aCij (3) 

h=£Aij • Yij (4) 

in which i=1, 2, ... t and j=1, 2, 
so the apparatus comprising: 

1) a storage means for storing information representing the interpolation coefficients Bij and Cij, 

2) an input means for inputting the coefficient a as one of real numbers including a range smaller than 0 and/or 
a range larger than 1 , the coefficient a determining the sharpness of a secondary image, which is reproduced 

55 from the interpolation image signal, 

3) an interpolation coefficient operation means for calculating the interpolation coefficient Aij in accordance 
with the coefficient a, the calculation being carried out from the interpolation coefficients Bij and Cij, which are 
stored in the storage means, and the coefficient a inputted from the input means, and 

4) an interpolating operation means for storing Formula (4) as the operation formula, and calculating the value 
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of an interpolated image signal component Y p , which corresponds to an interpolation point Xp. with Formula 
(4) in accordance with the interpolation coefficient Aij, which has been calculated by the interpolation coeffi- 
cient operation means, and each of the original image signal components Yij. 

5 In the interpolating operation apparatus for an image signal in accordance with the present invention, as for the two 

different interpolating functions for obtaining two interpolation images having different levels of sharpness, a B spline 
interpolating operation function should preferably be employed as the interpolating function corresponding to an image 
having a comparatively low sharpness, and a cubic spline interpolating operation function should preferably be 
employed as the interpolating function corresponding to an image having a sharpness higher than the comparatively 

w low sharpness. 

However, the interpolating operation apparatus for an image signal in accordance with the present invention is not 
limited to the combination of the B spline interpolating operation function and the cubic spline interpolating operation 
function. For example, various interpolating operation functions, such as the B spline interpolating operation function, 
the cubic spline interpolating operation function, a linear interpolating function, and a Lagrangean interpolating opera- 
75 tion function, may be used, and an arbitrary combination of two of these functions may be employed. 

The interpolating operation apparatus for an image signal in accordance with the present invention may be modi- 
fied such that the input means may be provided with: 

a) a response input means for inputting a response R desired for the interpolation image, 
20 b) a plurality of first look-up tables, each of the first look-up tables setting the frequencies and the corresponding 
responses with respect to one of the two interpolating functions and for each of a plurality of different image size 
enlargement scale factors, 

c) a plurality of second look-up tables, each of the second look-up tables setting the frequencies and the corre- 
sponding responses R 2 with respect to the other interpolating function and for each of a plurality of different image 

25 size enlargement scale factors, and 

d) a coefficient calculating means for calculating a response of the one interpolating function with respect to an 
image size enlargement scale factor, which is desired for an interpolation image, and a response R 2 of the other 
interpolating function with respect to the image size enlargement scale factor, which is desired for the interpolation 
image, by making reference to a first look-up table and a second look-up table, which correspond to the image size 

30 enlargement scale factor desired for the interpolation image, and 

for calculating the coefficient a from an operation carried out with Formula (5): 

«=(R-Ri)/(R 2 - R i) (5) 

35 in accordance with the response R desired for the interpolation image, the calculated response R-t of the one inter- 
polating function, and the calculated response R 2 of the other interpolating function. 

In such cases, the image size enlargement scale factor, which is the enlargement scale factor of the interpolation 
image with respect to the original image, may be inputted from an independent enlargement scale factor input means, 
40 or the like. Alternatively, the response input means may also serve as the means for inputting the image size enlarge- 
ment scale factor. 

As the first look-up table and the second look-up table corresponding to a predetermined image size enlargement 
scale factor, for example, those illustrated in Figures 5A and SB may be employed. As in the aforesaid interpolating 
operation method in accordance with the present invention, in cases where the first look-up table and the second look- 

45 up table corresponding to the inputted image size enlargement scale factor are not contained in the plurality of the look- 
up tables, two first look-up tables corresponding to two image size enlargement scale factors, which are closest to the 
inputted image size enlargement scale factor, may be selected from the first look-up table group. A response may then 
be calculated by carrying out an interpolating operation using an equation of the first order on the responses found from 
the two selected first look-up tables. Also, two second look-up tables corresponding to two image size enlargement 

so scale factors, which are closest to the inputted image size enlargement scale factor, may be selected from the second 
look-up table group. A response may then be calculated by carrying out an interpolating operation using an equation of 
the first order on the responses found from the two selected second look-up tables. 

The present invention further provides a first spline interpolating operation method for an image signal, wherein an 
interpolating operation using a third-order spline interpolating function is carried out on an original image signal, which 

55 represents an original image and is made up of a series of original image signal components, and an interpolation 
image signal is obtained with Formula (6), the interpolation image signal being made up of a series of image signal 
components, which occur at intervals different from those of the original image signal components, 

Y p=3 k-1 Y k-1 +a k Y k +a k+1 Y k+1 +a k+2 Y k+2 (®) 
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in which Y k _ 1( Y k , Y k+1 , and Y k+2 represent the original* image signal components representing picture elements X^, 
X k , and ^+2 in tne original image, Y p represents the interpolated image signal component corresponding to an 
interpolation point Xp located between the picture elements X k and X k+1 , and a^ , a k , a k+ i , and a k+ 2 represent the inter- 
polation coefficients, 

the method comprising the steps of: 

i) calculating the interpolation coefficients a k . 1( a^ a k+1 , and a k+2 respectively corresponding to the original image 
signal components Y k . v Y k , Y k+1 , and Y k+2 such that, in cases where the original image signal components repre- 
senting the two picture elements X k ^ and X k located before the interpolation point Xp, which is located between the 
picture elements X k and X k+1 in the original image, are represented by Y k .-i and Y k , the original image signal com- 
ponents representing the two picture elements X k+1 and X k+2 located after the interpolation point Xp are repre- 
sented by Y k+1 and Y k+2 , the third-order spline interpolating function between the picture elements X k and X k+1 is 
represented by f k , the third-order spline interpolating function between the picture elements X^ and X k is repre- 
sented by f k . 1t and the third-order spline interpolating function between the picture elements X k+1 and Xk +2 is rep- 
resented by f k+1 : 

a) the spline interpolating function f k at the picture elements X k and X k+1 may satisfy the original image signal 
components Y k and Y k+1 , as represented by Formulas (7) and (8): 

MX k )=Y k (7) 

M* k+ i)=Y k+1 (8) 

b) the first-order differential coefficient of the spline interpolating function f k at the picture element X k may coin- 
cide with the first-order differential coefficient of the spline interpolating function f k _! at the picture element X k , 
as represented by Formula (9): 

V(X k )=f k V(X k ) .(9) 

c) the first-order differential coefficient of the spline interpolating function f k at the picture element X k+1 may 
coincide with the first-order differential coefficient of the spline interpolating function f k+1 at the picture element 
X k+1 , as represented by Formula (10): 

^'(X k+1 )=f k+ i , (X k+1 ) (10) 

d) the first-order differential coefficient of the spline interpolating function f k at the picture element X k may have 
an inclination of an arbitrary parameter a with respect to the gradient of the original image signal components 
Y k ^ and Y k+1 representing the picture elements X^ and Xk +1 , which are located before and after the picture 
element X^ as represented by Formula (11), the arbitrary parameter a being selected previously and deter- 
mining the sharpness of a secondary image represented by the interpolation image signal, 

f k (X k )= a(Y k+1 -Y k ., )/(X k+1 -X k ., ) (11) 

and 

e) the first-order differential coefficient of the spline interpolating function f k at the picture element X k+1 may 
have an inclination of the parameter a with respect to the gradient of the original image signal components Y k 
and Y k+2 representing the picture elements X k and X^, which are located before and after the picture element 
X^, as represented by Formula (12): 

* k'(X k+1 )=a(Y k+2 -Y k )/(X k+2 -X k ) (12) 

and 

ii) calculating the interpolated image signal component Y p corresponding to the interpolation point Xp in accord- 
ance with the calculated interpolation coefficients a k . 1( a*, a k+1 , and a k+2 , and the original image signal compo- 
nents Y k . 1f Y k , Y k+1l and Y k+2 . 

Specifically, in the first spline interpolating operation method for an image signal in accordance with the present 
invention, the calculated interpolation coefficients a k . 1 , a*, a^, and a k+2 may respectively be represented by Formulas 
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(13). (14), (15), and (16): 

a K . 1= -(a/2)t(M) 2 (13) 

a k =(2-a/2)t 3 -(3-a/2)t 2 +1 (14) 

a k+1 =(o t /2-2)t 3 +(3-a)t 2 +(a/2)t (15) 

a k+2 =(a/2)t 2 (t-1) (16) 

in which t, where 0^t<1 . represents the position of the interpolation point Xp, the position being taken with respect to 
the picture element X^ that serves as a reference position, and in the direction heading toward the picture element 
Xk +1 , the lattice interval of the original image signal components (i.e., the picture element interval) being set to be equal 
to 1. 

The present invention still further provides a second spline interpolating operation method for an image signal, 
wherein an interpolating operation using a third-order spline interpolating function is carried out on an original image 
signal, which represents an original image and is made up of a series of original image signal components, and an inter- 
polation image signal is obtained with Formula (6), the interpolation image signal being made up of a series of image 
signal components, which occur at intervals different from those of the original image signal components, 

Y p =a k-1 Y k-1 +a k Y k +a k+1 Y k+1 +a k+2 ^ k+2 ( 6 ) 

in which Y k . 1( 'Y k , Y k+1 , and Y k+2 represent the original image signal components representing picture elements X k . 1t 
Xk. X k+1t and X^ in the original image, Y p represents the interpolated image signal component corresponding to an 
interpolation point Xp located between the picture elements X k and X k+1( and a^-, , a k , a^ , and a k+2 represent the inter- 
polation coefficients, 

the method comprising the steps of: 

i) calculating the interpolation coefficients a k . 1( ak, a^, and a^ respectively corresponding to the original image 
signal components Y k . 1f Y k , Y k+1 , and Y k+2 such that, in cases where the original image signal components repre- 
senting the two picture elements X^., and X k located before the interpolation point Xp, which is located between the 
picture elements X k and X k+1 in the original image, are represented by Y,^ and Y k , the original image signal com- 
ponents representing the two picture elements X^ and X k+2 located after the interpolation point Xp are repre- 
sented by Y k+1 and Y k+2l the third-order spline interpolating function between the picture elements Xk and X k+1 is 
represented by f k , the third-order spline interpolating function between the picture elements X k _-, and X k is repre- 
sented by f k . 1( and the third-order spline interpolating function between the picture elements X k+1 and X^ is rep- 
resented by f k+1 : 

a) the spline interpolating function f k at the picture elements X k and X k+1 may slightly deviate from the original 
image signal components Y k and Y k+1 in accordance with an arbitrary parameter p, as represented by Formu- 
las (1 7) and (18), the arbitrary parameter p being selected previously and determining the sharpness of a sec- 
ondary image represented by the interpolation image signal, 

f k(* k )=-0.5p Y m +0 + P)Y k -0.5p Y k+1 (17) 
* k(* k + i >=-0.5p Y k +(1 +P)Y k+1 -0.5PY k+2 (18) 

b) the first-order differential coefficient of the spline interpolating function f k at the picture element Xk may coin- 
cide with the first-order differential coefficient of the spline interpolating function f k ^ at the picture element Xk, 
as represented by Formula (9): 

WkHk-iW (9) 

c) the first-order differential coefficient of the spline interpolating function f k at the picture element Xk +1 may 
coincide with the first-order differential coefficient of the spline interpolating function f k+1 at the picture element 
Xk +1 , as represented by Formula (10): 

V(*k+i)=W(Xk + i) (10) 
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d) the first-order differential coefficient of the spline interpolating function f k at the picture element Xj< may coin- 
cide with the gradient of the original image signal components and Y k+1 representing the picture elements 

and X k+1 , which are located before and after the picture element X^ as represented by Formula (19): 

5 fkXXfcMY^-Y^WX^-X^) (19) 

and 

e) the first-order differential coefficient of the spline interpolating function f k at the picture element Xk +1 may 
coincide with the gradient of the original image signal components Y k and Y k+2 representing the picture ele- 

w ments and X k+2 , which are located before and after the picture element X k+1 , as represented by Formula 

(20): 

^'(X k+1 )=(Y k+2 -Y k )/(X k+2 -X k ) (20) 

is and 

ii) calculating the interpolated image signal component Y p corresponding to the interpolation point X p in accord- 
ance with the calculated interpolation coefficients a^, a k , a k+1 , and a k+2 , and the original image signal compo- 
nents Y k . 1( Y kl Y k+1 , and Y k+2 . 

20 

Specifically, in the second spline interpolating operation method for an image signal in accordance with the present 
invention, the calculated interpolation coefficients a k . 1( a k . a k+1 , and a^ may respectively be represented by Formulas 
(21), (22), (23). and (24): 

25 a k _ t =-(P+1fc)t 3 +(3p/2+1)t 2 -(1/2)t-p/2 (21) 

a k =3(p+1/2)t 3 -(9p/2+5/2)t 2 +p+1 (22) 
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a k+1 =-3(p+1/2)t 3 +(9p/2+2)t 2 +(1/2)t-p/2 (23) 
a k+2 =(p+1/2)t 3 -(3p/2+1/2)t 2 (24) 



in which t, where 0^t<1 , represents the position of the interpolation point Xp, the position being taken with respect to 
the picture element X k , that serves as a reference position, and in the direction heading toward the picture element 
35 Xk+1. the lattice interval of the original image signal components (i.e., the picture element interval) being set to be equal 
to 1. 

The present invention also provides a third spline interpolating operation method for an image signal, wherein an 
interpolating operation using a third-order spline interpolating function is carried out on an original image signal, which 
represents an original image and is made up of a series of original image signal components, and an interpolation 
40 image signal is obtained with Formula (6), the interpolation image signal being made up of a series of image signal 
components, which occur at intervals different from those of the original image signal components, 

^ p sa k-1 Y k-1 +a k +a k+1 Y k+1 +a k+2 Y k+2 ( 6 ) 

45 in which Y k . 1f Y k , Y k+1 , and Y k+2 represent the original image signal components representing picture elements X k . 1 , 
^k. and *k+2 in toe original image, Y p represents the interpolated image signal component corresponding to an 
interpolation point Xp located between the picture elements X k and X k+1 , and a k . 1 , ak, a k+1 , and a k+2 represent the inter- 
polation coeff id ents, 

the method comprising the steps of: 

50 

i) calculating the interpolation coefficients a^-,, a k , a k+1 . and a k+2 respectively corresponding to the original image 
signal components Y^ , Y k , Y k+1i and Y k+2 such that, in cases where the original image signal components repre- 
senting the two picture elements and Xj< located before the interpolation point Xp, which is located between the 
picture elements Xj< and X k+1 in the original image, are represented by Y k .-, and Y k , the original image signal com- 
55 ponents representing the two picture elements X k+1 and X k+2 located after the interpolation point Xp are repre- 
sented by Y k+1 and Y k+2 , the third-order spline interpolating function between the picture elements Xj< and Xk +1 is 
represented by f k , the third-order spline interpolating function between the picture elements X^ and X k is repre- 
sented by f k . 1t and the third-order spline interpolating function between the picture elements X|< +1 and X k+2 is rep- 
resented by f k+1 : 
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a) the spline interpolating function f k at the picture elements X k and X k+1 may satisfy the original image signal 
components Y k and Y k+1 , as represented by Formulas (7) and (8): 

MX k )=Y k (7) 
MX k+ i)=Y k+1 (8) 

b) the first-order differential coefficient of the spline interpolating function f k at the picture element X k may coin- 
cide with the first-order differential coefficient of the spline interpolating function f k-1 at the picture element X k , 
as represented by Formula (9): 

V(X k )=W(Xk) (9) 

c) the first-order differential coefficient of the spline interpolating function f k at the picture element X k+1 may 
coincide with the first-order differential coefficient of the spline interpolating function f k+1 at the picture element 
X k+1 , as represented by Formula (10): 

V(x k+1 )=f k+1 '(x k+l ) (10) 

d) the first-order differential coefficient of the spline interpolating function f k at the picture element may have 
an inclination of an arbitrary parameter a with respect to the gradient of the original image signal components 
Y k _i and Y k+1 representing the picture elements X k . t and X k+1 , which are located before and after the picture 
element X k , as represented by Formula (11), the arbitrary parameter a being selected previously and deter- 
mining the sharpness of a secondary image represented by the interpolation image signal, 

f k '(X k )=a(Y k+1 -Y k ., )/(X k+1 -X k-1 ) (11) 

and 

e) the first-order differential coefficient of the spline interpolating function f k at the picture element X k+1 may 
have an inclination of the parameter a with respect to the gradient of the original image signal components Y k 
and Y k+2 representing the picture elements X k and X k+2 , which are located before and after the picture element 
X k+1 , as represented by Formula (12): 

f k '(X k+1 )=ct(Y k+2 -Y k )/(X k+2 -X k ) (12) 

ii) calculating interpolation coefficients b^, t^, b k+1 , and b^ +2 respectively corresponding to the original image sig- 
nal components Y k-1l Y k , Y k+1 , and Y k+2 such that: 

0 the spline interpolating function f k at the picture elements and X^-, may slightly deviate from the original 
image signal components Y k and Y k+1 in accordance with an arbitrary parameter p, as represented by Formu- 
las (17) and (18), the arbitrary parameter p being selected previously and determining the sharpness of the 
secondary image represented by the interpolation image signal, 

f k (X k )=-0.5pY k . 1+ (1 + p)Y k -0.5pY k+1 (17) 

* k(* k + i )~0.5P Y k+ (1 +P)Y k+1 -0.5PY k+2 (1 8) 

g) the first-order differential coefficient of the spline interpolating function f k at the picture element X k may coin- 
cide with the first-order differential coefficient of the spline interpolating function f k-1 at the picture element X^ 
as represented by Formula (9): 

V<x k M*r<x k ) (9) 

h) the first-order differentia! coefficient of the spline interpolating function f k at the picture element X k+1 may 
coincide with the first-order differential coefficient of the spline interpolating function f k+1 at the picture element 
X k+1 , as represented by Formula (10): 

V(X k+l )=W(X k+1 ) (10) 
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i) the first-order differential coefficient of the spline interpolating function f k at the picture element may coin- 
cide with the gradient of the original image signal components and Y k+1 representing the picture elements 
and Xk+v which are located before and after the picture element X|<, as represented by Formula (19): 

5 ^XX^CYk^-Y^JAX^-X^) (19) 

and 

j) the first-order differential coefficient of the spline interpolating function f k at the picture element X|< +1 may 
coincide with the gradient of the original image signal components Y k and Y k+2 representing the picture ele- 
w ments X k and Xk+2. which are located before and after the picture element X k+ i . as represented by Formula 

(20): 

^ , (X k+1 )=(Y k+r Y k )/(X k+2 -X k ) (20) 

is iii) calculating mean values (a ^ +b ^ )/2 , (a k +b k )/2 , (a k+1 +b k+1 )/2 , and (a k+2 +b k+2 )/2 of the respective sets 
of the interpolation coefficients, which sets respectively correspond to the original image signal components Y k . 1( 
Y k , Y k+1 , and Y k+2 . the calculated mean values being taken as the values of new interpolation coefficients a k .j, a k , 
a k+1 , and a k+2 , and 

iv) calculating the interpolated image signal component Y p corresponding to the interpolation point X p in accord- 
20 ance with the calculated new interpolation coefficients a^i, ak. a k+1t and ak+ 2 , and the original image signal com- 
ponents Yk.! , Y k , Y k+1 , and Y k+2 . 

Specifically, in the third spline interpolating operation method for an image signal in accordance with the present 
invention, the calculated new interpolation coefficients a k . 1( a k , ak +1 , and ak +2 , which represent the mean values, may 
25 respectively be represented by Formulas (25), (26), (27), and (28): 

a k _ 1 H(a+2p+1)/4}t 3 +{(2a+3p+2)/4}t 2 -{(a+1)/4}t-p/4 (25) 



30 



a k ={(-a+6p+7)/4}t 3 +{(a-9p-11)/4}t 2 +P/2+1 (26) 

a k+1 ={(a-6p-7)/4}t 3 +{(-2a+9p+10)/4}t 2 +{(a+1)/4}t-p/4 (27) 

a k+2 = {( a +2 p + 1 )/4}t 3 -{( a+3 p + 1 )/4}t 2 (28) 

35 in which t, where 0^t<1 , represents the position of the interpolation point Xp, the position being taken with respect to 
the picture element Xk, that serves as a reference position, and in the direction heading toward the picture element 
Xk + i, the lattice interval of the original image signal components (i.e., the picture element interval) being set to be equal 
to 1. 

The present invention further provides a first spline interpolating operation apparatus for an image signal, wherein 
40 an interpolating operation using a third-order spline interpolating function is carried out on an original image signal, 
which represents an original image and is made up of a series of original image signal components, and an interpola- 
tion image signal is obtained with Formula (6), the interpolation image signal being made up of a series of image signal 
components, which occur at intervals different from those of the original image signal components, 

45 Y p =a k _., Y k _.,+a k Y k +a k+1 Y k+1 +a k+2 Y k+2 (6) 

in which Y k _ 1( Y k , Y k+1 , and Y k+2 represent the original image signal components representing picture elements X^, 
Xk, Xk+i , and X|< +2 in the original image, Y p represents the interpolated image signal component corresponding to an 
interpolation point Xp located between the picture elements Xk and Xk+i. and ak-i. ak, a k+ i. and ak+2 represent the inter- 
so polation coefficients, 

the apparatus comprising: 

1) a storage means for storing information representing the interpolation coefficients a^, a k , a^-i. and ak+2. which 
respectively correspond to the original image signal components Y k . 1( Y k , Y k+1 , and Y k+2 representing the picture 
55 elements X^, X k , Xk+i, and Xk+ 2 in the original image, and which are set to be respectively represented by For- 
mulas (13). (14), (15), and (16): 

a K . 1= -(a/2)t(t-1) 2 (13) 
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a k =(2-a/2)t 3 -(3-a/2)t?+1 (14) 

a k+1 =(o/2-2)t 3 +(3-a)t 2 +(a/2)t { 1 5) 

5 a k+2 =(a/2)t 2 (t-1) (16) 

in which a represents the parameter determining the sharpness of a secondary image represented by the interpo- 
lation image signal, and t. where 0^t<1, represents the position of the interpolation point X p , the position being 
taken with respect to the picture element X^ that serves as a reference position, and in the direction heading 
10 toward the picture element X k+1( the lattice interval of the original image signal components being set to be equal 
to1. 

2) an input means for inputting the arbitrary parameter a, which determines the sharpness of the secondary image 
reproduced from the interpolation image signal, 

3) an interpolation coefficient operation means for calculating the interpolation coefficients a^-,, a k , a k+1f and a k+2 
is in accordance with the parameter a, the calculation being carried out from the interpolation coefficients, which are 

stored in the storage means, and the parameter a inputted from the input means, and 

4) an interpolating operation means for storing Formula (6) as the operation formula, and calculating the value of 
the interpolated image signal component Y pj which corresponds to the interpolation point X p , with Formula (6) in 
accordance with the interpolation coefficients a k ^ , a k , a k+1 , and a k+2 . which have been calculated by the interpola- 
te tion coefficient operation means, and the original image signal components Y k . 1f Y k , Y k+1 , and Y k+2 . 

The present invention still further provides a second spline interpolating operation apparatus for an image signal, 
wherein an interpolating operation using a third-order spline interpolating function is carried out on an original image 
signal, which represents an original image and is made up of a series of original image signal components, and an inter- 
ns polation image signal is obtained with Formula (6), the interpolation image signal being made up of a series of image 
signal components, which occur at intervals different from those of the original image signal components, 

Y p =a k-1 Y k-1 +a k Y k +a k+i Y k+1 +a k+2 Y k+2 ( 6 ) 

30 in which Y k . 1p Y k , Y k+1 , and Y k+2 represent the original image signal components representing picture elements X^, 
X^ X^, and X k+2 in the original image, Y p represents the interpolated image signal component corresponding to an 
interpolation point X p located between the picture elements and Xk +1 , and a k . 1 , a*, a^ , and a k+2 represent the inter- 
polation coefficients, 

the apparatus comprising: 

35 

1) a storage means for storing information representing the interpolation coefficients a^, a k , a k+1 , and a k+2 , which 
respectively correspond to the original image signal components Y k . 1f Y k , Y k+1 , and Y k+2 representing the picture 
elements X k . 1 , X k , X k+1 , and Y^+ 2 in tne original image, and which are set to be respectively represented by For- 
mulas (21), (22), (23), and (24): 

40 

a k-i=-(P+' | / 2 )t 3 +(3p/2+1)t 2 -(1/2)t-p/2 (21) 

a k =3(p+1/2)t 3 -(9p/2+5/2)t 2 +p+1 (22) 

« a k+1 =-3(p+1/2)t 3 +(9p/2+2)t 2 +(1/2)t-p/2 (23) 

a k+2 =(p+1/2)t 3 -(3p/2+1/2)t 2 (24) 

in which p represents the parameter determining the sharpness of a secondary image represented by the interpo- 
se? lation image signal, and t, where 0^t<1 , represents the position of the interpolation point Xp, the position being 
taken with respect to the picture element X^ that serves as a reference position, and in the direction heading 
toward the picture element X k+1 , the lattice interval of the original image signal components being set to be equal 
to 1, 

2) an input means for inputting the arbitrary parameter p, which determines the sharpness of the secondary image 
55 reproduced from the interpolation image signal, 

3) an interpolation coefficient operation means for calculating the interpolation coefficients a^. a^ a k+1 . and a k+2 
in accordance with the parameter p, the calculation being carried out from the interpolation coefficients, which are 
stored in the storage means, and the parameter p inputted from the input means, and 

4) an interpolating operation means for storing Formula (6) as the operation formula, and calculating the value of 
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the interpolated image signal component Y p , which .corresp9nds to the interpolation point Xp, with Formula (6) in 
accordance with the interpolation coefficients a k . 1( a k , a k+1 , and a^, which have been calculated by the interpola- 
tion coefficient operation means, and the original image signal components Y k _ 1( Y k , Y k+1 , and Y k+2 . 

5 The present invention also provides a third spline interpolating operation apparatus for an image signal, wherein an 

interpolating operation using a third-order spline interpolating function is carried out on an original image signal, which 
represents an original image and is made up of a series of original image signal components, and an interpolation 
image signal is obtained with Formula (6), the interpolation image signal being made up of a series of image signal 
components, which occur at intervals different from those of the original image signal components, 

10 

Y p =a k-1 Y k-1 +a k Y k +a k+1 Y k+1 +a k+2 Y k+2 ( 6 ) 

in which Y k _ 1( Y k( Y k+1 , and Y k+2 represent the original image signal components representing picture elements X k . 1t 
x k- Xk+L and x k+2 in the original image, Y p represents the interpolated image signal component corresponding to an 
15 interpolation point Xp located between the picture elements Xj< and X k+1 , and a k . 1 , a k , a k+1 , and a k+2 represent the inter- 
polation coefficients, 

the apparatus comprising: 

1) a storage means for storing information representing the interpolation coefficients a^-j, a k , a k+1 . and a k+2 , which 
20 respectively correspond to the original image signal components Y k . v Y k , Y k+1 , and Y k+2 representing the picture 
elements X^-,, X^ X k+1 , and X|< +2 in the original image, and which are set to be respectively represented by For- 



mulas (25), (26), (27), and (28): 

a k . 1 =-{(a+2P+1)/4}t 3 +{(2a+3p+2)/4}t 2 -{(a+1)/4}t-p/4 (25) 

25 

a K ={(-a+6p+7)/4}t 3 +{(a-9p-11)/4}t 2 +P/2+1 (26) 

a k+1 ={(a-6p-7)/4}t 3 +{(-2a+9p+10)/4}t 2 +{(a+1)/4}t-p/4 (27) 

30 a k+2 ={(a+2p+1)/4}t 3 -{(a+3P+1)/4}t 2 (28) 



in which a and p represent the parameters determining the sharpness of a secondary image represented by the 
interpolation image signal, and t, where 0^t<1. represents the position of the interpolation point Xp, the position 
being taken with respect to the picture element X kl that serves as a reference position, and in the direction heading 
35 toward the picture element X^, the lattice interval of the original image signal components being set to be equal 
to 1. 

2) an input means for inputting the arbitrary parameters a and p, which determine the sharpness of the secondary 
image reproduced from the interpolation image signal, 

3) an interpolation coefficient operation means for calculating the interpolation coefficients a^. a^ a k+1 , and a k+2 
40 in accordance with the parameters a and p, the calculation being carried out from the interpolation coefficients, 

which are stored in the storage means, and the parameters a and p inputted from the input means, and 

4) an interpolating operation means for storing Formula (6) as the operation formula, and calculating the value of 
the interpolated image signal component Y p , which corresponds to the interpolation point Xp, with Formula (6) in 
accordance with the interpolation coefficients a k . n , a k , a^, and a k+2 , which have been calculated by the interpola- 

45 tion coefficient operation means, and the original image signal components Y k . 1( Y k , Y k+1 , and Y k+2 . 

With the interpolating operation method and apparatus for an image signal in accordance with the present inven- 
tion, the new interpolation coefficient is obtained for each of the original image signal components of the original image 
signal representing the original image by linearly combining the interpolation coefficients, which correspond to each 

so other and are set for each of the original image signal components, in a first interpolating function f for obtaining an 
interpolation image having a first sharpness and a second interpolating function g for obtaining an interpolation image 
having a second sharpness different from the first sharpness. At this time, the coefficient for the weighting of the two 
interpolation coefficients in the linear combination is not limited to the range of 0 to 1 . and is set to be an arbitrary real 
number. In this manner, interpolation images having a wide variety of levels of sharpness, which are not limited to the 

55 sharpness range between the first sharpness and the second sharpness, can be obtained. 

Specifically, the first interpolating function f and the second interpolating function g may be represented by the for- 
mulas 

f=£Bij-Yij, g=£Cij-Yij 
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in which Yij represents the original image signal components, and Bij and Cij represent the interpolation coefficients. In 
such cases, the term "linearly combining interpolation coefficients, which correspond to each other and are set for each 
of original image signal components" as used herein means that the operation is carried out with, for example, the for- 
mula 

A12=(1-ct)B12+aCl2 

wherein a represents all real numbers, on the interpolation coefficient B12 in the interpolating function f and the inter- 
polation coefficient C12 in the interpolating function g. which coefficients are for the original image signal component 
Y12 (in cases where i=1, and j=2). 

Also, as for the interpolation coefficients Bij and Cij for the other original image signal components Yij ( the interpo- 
lation coefficients Bij and Cij are linearly combined with each other by carrying out the operations with Formula (3) 

Aij=(1-a)Bij+aCij (3) 

and the new interpolation coefficients Aij are thereby obtained. 

In the conventional technique for the linear combination, the coefficient a was limited to the real numbers falling 
within the range of 0 to 1 . However, with the interpolating operation method and apparatus for an image signal in 
accordance with the present invention, the coefficient a can take even a value less than 0 and a value more than 1. 
Therefore, an interpolation image can be obtained, which has a sharpness higher than the sharpness of the interpola- 
tion image obtained with the interpolating function f or g ( whichever is capable of yielding an interpolation image having 
a higher sharpness. Alternatively, an interpolation image can be obtained, which has a lower sharpness and is 
smoother than the interpolation image obtained with the interpolating function f or g, whichever is capable of yielding 
an interpolation image, that has a lower sharpness and is smoother. Accordingly, the degree of freedom of the selection 
of sharpness of the interpolation image in accordance with the kind of the image and the image size enlargement scale 
factor can be kept markedly high. For example, as for a radiation image of a blood vessel, a very sharp image is desired. 
In such cases, the coefficient for the weighting of the interpolation coefficient (e.g.. Cij) in the interpolating function, 
which yields a higher sharpness, may be set to be a value larger than 1. As a result, an interpolation image having a 
very high sharpness can be obtained. As for a CT scanner image of the liver, it is desired that a step-wise density 
change portion, which often occur due to a low CT resolution, can be reproduced to be unsharp. In such cases, the 
coefficient for the weighting of the interpolation coefficient (e.g., Bij) in the interpolating function, which yields a lower 
sharpness, may be set to be a value larger than 1 . As a result, an interpolation image having a very low sharpness can 
be obtained. In this manner, the interpolating operation method and apparatus for an image signal in accordance with 
the present invention can meet various demands. 

In the interpolating operation apparatus for an image signal in accordance with the present invention, the input 
means may be provided with the response input means for inputting a response R desired for the interpolation image, 
a plurality of the first look-up tables with respect to one of the two interpolating functions, a plurality of the second look- 
up tables with respect to the other interpolating function, and the coefficient calculating means. In such cases, the 
response R of the image, with which the extent of the change in sharpness can be grasped easily, can be specified 
directly as a desired value. The coefficient a corresponding to the response is then calculated by using the first and sec- 
ond look-up tables. Therefore, an interpolation image having a sharpness truer to nature can be obtained than with the 
apparatus for specifying the value of a serving as a mere inorganic coefficient. 

With the first spline interpolating operation method and apparatus for an image signal in accordance with the 
present invention, the first-order differential coefficient of the third-order spline interpolating function f k corresponding to 
the section X^X^-,, in which the interpolation point Xp is located, at the picture element ^ is set to have an inclination 
of an arbitrary parameter a with respect to the gradient of the original image signal components and Y k+1 repre- 
senting the picture elements X^ and X^, which are located before and after the picture element Xk, the arbitrary 
parameter a being selected previously and determining the sharpness of a secondary image represented by the inter- 
polation image signal. Also, the first-order differential coefficient of the spline interpolating function f k at the picture ele- 
ment X k+1 is set to have an inclination of the parameter a with respect to the gradient of the original image signal 
components Y k and Y k+2 representing the picture elements X,< and X^, which are located before and after the picture 
element X|< +1 . Therefore, an interpolation image signal having a desired level of sharpness can be obtained by arbitrar- 
ily altering the value of the parameter a. 

Specifically, with the conventional cubic spline interpolating operation, as represented by Formulas (19) and (20), 
the parameter a in Formulas (11) and (12) is fixed to be equal to 1. Therefore, in cases where the conventional cubic 
spline interpolating operation is carried alone, or in cases where the conventional B spline interpolating operation is car- 
ried out alone, only the secondary image having a single level of sharpness can be obtained. Also, with the spline inter- 
polating operation, in which the cubic spline interpolating operation and the B spline interpolating operation are 
combined with each other and which is disclosed in U.S.R 5.048,105, the adjustment of sharpness can be carried out 
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only within the range of the sharpness, which is obtained .with thejcubic spline interpolating operation, to the sharpness, 
which is obtained with the B spline interpolating operation. However, with the first spline interpolating operation method 
for an image signal in accordance with the present invention, by the setting of the parameter a to be larger than 1 , an 
interpolation image having a sharpness higher than the sharpness of the interpolation image, which is obtained from 
the conventional cubic spline interpolating operation, can be obtained. As the parameter a is set to be a larger value, a 
higher sharpness can be obtained. Also, by the setting of the parameter a to be smaller than 1 . an interpolation image 
having a sharpness lower than the sharpness of the interpolation image, which is obtained from the conventional B 
spline interpolating operation, can be obtained. 

As the parameter a is set to be a smaller value, a lower sharpness can be obtained. 

Further, the calculation processing with the first spline interpolating operation method for an image signal in accord- 
ance with the present invention can be carried out only with the same calculation processing as the conventional cubic 
spline interpolating operation. Therefore, with the first spline interpolating operation method for an image signal in 
accordance with the present invention, the calculation processing can be kept simpler than in the spline interpolating 
operation disclosed in U.S.R 5,048,105. Accordingly, the time required to carry out the calculation processing can be 
kept short, while the adjustment of sharpness is being enabled. 

With the second spline interpolating operation method and apparatus for an image signal in accordance with the 
present invention, the values of the third-order spline interpolating function f k corresponding to the section X^X^ , in 
which the interpolation point Xp is located, at the picture elements X k and X^ are set to slightly deviate from the orig- 
inal image signal components Y k and Y k+1 in accordance with an arbitrary parameter p, which is selected previously 
and determines the sharpness of a secondary image represented by the interpolation image signal. Therefore, an inter- 
polation image signal having a desired level of sharpness can be obtained by arbitrarily altering the value of the param- 
eter p. 

Specifically, in cases where the conventional cubic spline interpolating operation is carried alone, or in cases where 
the conventional B spline interpolating operation is carried out alone, only the secondary image having a single level of 
sharpness can be obtained. Also, with the spline interpolating operation, in which the cubic spline interpolating opera- 
tion and the B spline interpolating operation are combined with each other and which is disclosed in U.S.R 5,048,105, 
the adjustment of sharpness can be carried out only within the range of the sharpness, which is obtained with the cubic 
spline interpolating operation, to the sharpness, which is obtained with the B spline interpolating operation. However, 
with the second spline interpolating operation method for an image signal in accordance with the present invention, by 
the setting of the parameter p to be equal to a value other than 0, an interpolation image having a sharpness higher 
than the sharpness of the interpolation image, which is obtained from the conventional cubic spline interpolating oper- 
ation, can be obtained. Also, an interpolation image having a sharpness lower than the sharpness of the interpolation 
image, which is obtained from the conventional B spline interpolating operation, can be obtained. 

Further, as represented by Formulas (1 7) and (1 8), in cases where the value on one end side of the spline interpo- 
lating function f k corresponding to the section X k ~X k+1 is increased, the value on the other end side of the spline inter- 
polating function f k is decreased. Therefore, even if the interpolating operation is carried out on an original image, the 
original image signal components of which are distributed with approximately equal values, the values of the interpola- 
tion image signal will not become far different from the values of the original image signal. Accordingly, an interpolation 
image having a natural feeling can be obtained. 

Furthermore, the calculation processing with the second spline interpolating operation method for an image signal 
in accordance with the present invention can be carried out only with the same calculation processing as the conven- 
tional cubic spline interpolating operation. Therefore, with the second spline interpolating operation method for an 
image signal in accordance with the present invention, the time required to carry out the calculation processing can be 
kept short, while the adjustment of sharpness is being enabled. 

With the third spline interpolating operation method and apparatus for an image signal in accordance with the 
present invention, the mean values of the values of the interpolation coefficients, which are obtained with the aforesaid 
first spline interpolating operation method (and apparatus) for an image signal, and the values of the corresponding 
interpolation coefficients, which are obtained with the aforesaid second spline interpolating operation method (and 
apparatus) for an image signal, are taken as the values of the new interpolation coefficients. The two spline interpolating 
operation methods (and apparatuses) can thus be combined with each other, and the degree of freedom of the sharp- 
ness adjustment can thereby be kept higher. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a schematic block diagram showing an image reproducing system provided with an embodiment of the 
interpolating operation apparatus for carrying out the interpolating operation method for an image signal in accord- 
ance with the present invention, 

Figure 2 is a perspective view showing an image read-out apparatus, 

Figure 3 is a block diagram showing an input means, which is provided with a response input means, first and sec- 
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ond look-up tables, and a coefficient calculating means, 

Figure 4 is an explanatory graph showing how interpolated image signal components are obtained with a conven- 
tional cubic spline interpolating operation from original image signal components, which are sampled with a period 
of an equal interval and represent sampling points (picture elements) arrayed in one direction, 
Figure 5A is a schematic graph showing an example of a first look-up table, the graph serving as an aid in explain- 
ing the relationship between a frequency and a response, 

Figure 5B is a schematic graph showing an example of a second look-up table, the graph serving as an aid in 
explaining the relationship between a frequency and a response, 

Figure 6 is a schematic block diagram showing an image reproducing system provided with an embodiment of the 
spline interpolating operation apparatus for carrying out the spline interpolating operation method for an image sig- 
nal in accordance with the present invention, 

Figure 7 is an explanatory graph showing how interpolated image signal components are obtained from original 
image signal components, which are sampled with a period of an equal interval and represent sampling points (pic- 
ture elements) arrayed in one direction, in the embodiment of the spline interpolating operation apparatus, 
Figures 8A, 8B, and 8C are graphs showing interpolation coefficients obtained when the value of a parameter a is 
set to be equal to 0.4, 1.0, and 1.6, respectively, 

Figures 9A, 9B ( and 9C are graphs showing interpolation coefficients obtained when the value of a parameter p is 
set to be equal to -0.5, 0. and 0.5, respectively, and 

Figures 10A, 10B, and 10C are graphs showing interpolation coefficients obtained when the values of the param- 
eters (a, p) are set to be equal to (0.4, -0.5), (1.0, 0.0), and (1.6, 0.5), respectively. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

The present invention will hereinbelow be described in further detail with reference to the accompanying drawings. 

Figure 1 is a schematic block diagram showing an image reproducing system provided with an interpolating oper- 
ation apparatus 30, which is an embodiment of the interpolating operation apparatus for carrying out the interpolating 
operation method for an image signal in accordance with the present invention. With reference to Figure 1, the image 
reproducing system comprises an image signal storing device 10 for storing an image signal Sorg representing an 
image, and a multi -formatter 20 for carrying out predetermined signal processing on the image signal (hereinbelow 
referred to as the primary image signal or the original image signal) Sorg, which is received from the image signal stor- 
ing device 10. such that an image signal conforming to a predetermined reproduction format may be obtained. The 
image reproducing system also comprises an image reproducing means 40 for reproducing a visible image in the pre- 
determined reproduction format from an image signal (hereinbelow referred to as the secondary image signal or the 
interpolation image signal) S\ which has been obtained from the predetermined signal processing carried out by the 
multi-formatter 20. 

The multi-formatter 20 carries out the signal processing on the primary image signal Sorg, such that an image sig- 
nal conforming to one of various image reproduction formats may be obtained. Examples of the image reproduction for- 
mats include a format, in which the entire area of a single sheet of film is divided into four different small regions and 
four different images having reduced image sizes are printed respectively in the four regions, a format, in which a single 
large image is printed directly on a single sheet of film, and a format, in which a portion of an image is enlarged and the 
enlarged image portion is printed on a sheet of film. The multi -formatter 20 is provided with the interpolating operation 
apparatus 30 in accordance with the present invention, which carries out an interpolating operation for obtaining the 
secondary image signal (the interpolation image signal) made up of a number of image signal components different 
from that of the image signal components of the primary image signal Sorg when the image size is to be enlarged or 
reduced. 

The primary image signal Sorg used in this embodiment is made up of a series of digital image signal components 
Yr-2. Y k . 1t Y kl Y k+1 , Y k+2 , ... respectively corresponding to sampling points (picture elements) X^. X k . 1( Xk. Xk +1 , Xk +2 , 
.... which are sampled with a period of an equal interval and arrayed in one direction. 

The interpolating operation apparatus 30, which is incorporated in the multi-formatter 20, comprises a cubic spline 
interpolation coefficient storing means 32, a B spline interpolation coefficient storing means 31, an interpolation coeffi- 
cient operation means 33, an input means 35. and an interpolating operation means 34. 

The cubic spline interpolation coefficient storing means 32 stores information representing interpolation coeffi- 
cients , Ck, C|< +1 . and c^. which respectively correspond to the original image signal components Y k _-, , Y k , Y k+1 , and 
Y k+2. in Formula (39). Formula (39) serves as a third-order cubic spline interpolating operation formula and.represents 
a first secondary image signal component Y p 1 corresponding to an interpolation point Xp, which is located between 
original sampling points Xk and X^. The interpolation coefficients c^i. c^ c k+1 , and q< + 2 are stored as being the ones 
represented by the formulas shown below. 

Y p 1 =c k-1 Y k-1 +° k Y k+1 Y k+1 + 0 k+2 Y k+2 ( 39 ) 
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C k . 1 =(-t?+2t 2 -t)^ t 

c k =(3t 3 -5t 2 +2)/2 

5 C k+1 =(-3t 3 +4t 2 +t)/2 

C k+2 =(t 3 -t 2 )/2 

in which t, where Ost£1, represents the position of the interpolation point Xp, the position being taken with respect to 
10 the picture element X k , that serves as a reference position, and in the direction heading toward the picture element 
X k+1 , the lattice interval being set to be equal to 1 . 

The B spline interpolation coefficient storing means 31 stores information representing interpolation coefficients b k . 
1. b k. bk+L and b k+2. which respectively correspond to the original image signal components Y k . 1f Y kl Y k+1 , and Y k+2 , 
in Formula (40). Formula (40) serves as a third-order B spline interpolating operation formula and represents a second 
15 secondary image signal component Y p 2 corresponding to the interpolation point Xp, which is located between original 
sampling points X k and X k+1 . The interpolation coefficients b k . 1( b k , b k+1 . and b k+2 are stored as being the ones repre- 
sented by the formulas shown below. 
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Y p 2=b M Y M +b k Y k+ b k+1 Y k+1 +b k+2 Y k+2 (40) 
b k . lB (-t 3 ^t 2 -3t+1)/6 
b k =(3t 3 -6t 2 +4)/6 

25 b k+1 =(-3t 3 +3t 2 +3t+1)/6 

b k+2 =t 3 /6 

in which t, where 0=*t==1. represents the position of the interpolation point Xp, the position being taken with respect to 
30 the picture element X k , that serves as a reference position, and in the direction heading toward the picture element 
X k+1 , the lattice interval being set to be equal to 1 . 

The interpolation coefficient operation means 33 receives the interpolation coefficients (hereinbelow referred to as 
the cubic spline interpolation coefficients) c^, c^, Ck +1 , and c k+2 from the cubic spline interpolation coefficient storing 
means 32 and receives the interpolation coefficients (hereinbelow referred to as the B spline interpolation coefficients) 
35 b k .t , b^ b k+1 , and bj< +2 from the B spline interpolation coefficient storing means 31 . The interpolation coefficient opera- 
tion means 33 weights and adds the cubic spline interpolation coefficients c^ . c k , q< +1 , and q< +2 and the B spline inter- 
polation coefficients b k . 1( b^. bk +1 , and b k+2l which coefficients correspond to each other and are set for the respective 
original image signal components Y k . 1( Y kl Y k+1 , and Y k+2 . The weighted additions are carried out with Formulas (41), 
(42), (43), and (44): 
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a^O-cOc^+abk.! (41) 



={(2a-3)t 3 -(3a-6)t 2 -3t+a}/6 



45 



55 



a k =(1-a)c k +ab k (42) 



={(9-6a)t 3 +(9a-15)t 2 +(6-2a)}/6 



a k+ i=( 1 - a ) c k + l +ab k + i ( 43 > 



50 9 

={(6a-9)t -(9a-12)t +3t+a}/6 

3 k + 2=( 1 " a ) c k + 2 +ab k+2 t 44 ) 
={(3-2a)t 3 +(3a-3)t 2 }/6 

The input means 35 is used to input an arbitrary value of the parameter a, which determines the weighting ratio, 
into the interpolation coefficient operation means 33. 

The interpolating operation means 34 stores Formula (38) shown below as the operation formula with the third- 
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order spline interpolating function and calculates the value of the interpolated image signal component Y p , which cor- 
responds to the interpolation point X p , with Formula (38) in accordance with the interpolation coefficients a^ , a*, a|< +1 , 
and afc +2 according to the parameter a, which have been calculated by the interpolation coefficient operation means 33, 
and the original image signal components Y k . 1f Y k , Y k+1 , and Y k+2 . 

^ p" a k-1 Y k-1 +a k Y k +a k+1 Y k+1 +a k+2^ k+2 (38) 

The cubic spline interpolation coefficients c^, c^, and c^, which are stored in the cubic spline interpolation 
coefficient storing means 32, and the B spline interpolation coefficients bk. 1( b k , b^, and b^, which are stored in the 
B spline interpolation coefficient storing means 31 , are calculated previously with the algorithms described above. Also, 
the parameter a can take all real numbers, including the range more than 1 and the range less than 0. 

An actual image is composed of the picture elements arrayed in two-dimensional directions. Therefore, the interpo- 
lation coefficients a k . 1t a^ a k+1l and a k+2 are calculated for each of two different directions (an i direction and a j direc- 
tion) of an array of the picture elements in the image. The interpolation coefficients a k _., , a k , a k+1 , and a k+2 having been 
calculated in this manner are often represented as the interpolation coefficients Aij. The B spline interpolation coeffi- 
cients b k . 1t b k , b k+1 , and b k+2 , which have been calculated for each of the i direction and the j direction of the array of 
the picture elements in the image, are often represented as the interpolation coefficients Bij. Also, the cubic spline inter- 
polation coefficients q^, c^ C|< +1 , and c k+2 , which have been calculated for each of the i direction and the j direction of 
the array of the picture elements in the image, are often represented as the interpolation coefficients Cij. 

In cases where the interpolation coefficient operation means 33 is caused to previously store Formulas (41), (42), 
(43), and (44), it can also serve as the B spline interpolation coefficient storing means 31 and the cubic spline interpo- 
lation coefficient storing means 32. 

In the image reproducing system provided with this embodiment, the interpolation image signal S' obtained from 
the interpolating operation means 34 may be directly fed into the image reproducing means 40. Alternatively, the inter- 
val of the image signal components of the interpolation image signal S* may be extended such that it may become equal 
to the interval of the image signal components of the original image signal Sorg. In this manner, the interpolation image 
can be reproduced as an image enlarged from the original image. The enlargement processing is carried out with the 
functions of the multi-formatter 20. For this purpose, the multi -formatter 20 receives information, which represents a 
desired image size enlargement scale factor, from an input means (not shown). 

How the image reproducing system provided with this embodiment operates will be described hereinbelow. 

Firstly, the multi-formatter 20 reads the primary image signal Sorg from the image signal storing device 10. Also, in 
order to obtain the secondary image signal, which represents an enlarged image corresponding to the image size 
enlargement scale factor inputted from the input means (not shown), the multi-formatter 20 feeds the primary image sig- 
nal Sorg into the interpolating operation apparatus 30 incorporated in the multi-formatter 20. 

In the interpolating operation apparatus 30, the primary image signal Sorg is fed into the interpolating operation 
means 34. 

The B spline interpolation coefficient storing means 31 and the cubic spline interpolation coefficient storing means 
32 set the values of t in the respective interpolation coefficients in accordance with the image size enlargement scale 
factor having been inputted from the input means (not shown), which is connected to the multi-formatter 20. For exam- 
ple, in cases where an image size enlargement scale factor of 2 is inputted, values of 0.5 and 1 .0 are set as the values 
of t In cases where an image size enlargement scale factor of 4 is inputted, values of 0.25, 0.5, 0.75, and 1.0 are set 
as the values of t. In cases where an image size enlargement scale factor of 10 is inputted, values of 0.1, 0.2, ... ( 1 .0 
are set as the values of t. Information representing the B spline interpolation coefficients and the cubic spline interpo- 
lation coefficients, which are for each of the thus set values of t f is fed into the interpolation coefficient operation means 
33. 

Also, information representing a value of the parameter (i.e., the coefficient) a corresponding to a level of sharp- 
ness desired for the secondary image is inputted into the input means 35. The information representing the value of the 
parameter a is also fed into the interpolation coefficient operation means 33. 

The information representing the value of the parameter a may be directly inputted by an operator from the exterior. 
Alternatively, information representing the response R. which corresponds to the level of sharpness desired for the 
interpolation image, may be inputted by the operator, and the response R may be converted into the parameter a in the 
input means 35. 

In order for this input means 35 to be provided with the functions for receiving the response R and converting it into 
the parameter a, the input means 35 may be provided with the constitution illustrated in Figure 3. 

Specifically, the input means 35 illustrated in Figure 3 is provided with a response input means 35a. from which the 
response R desired for the interpolation image is inputted by the operator. The input means 35 is also provided with six 
kinds of first look-up tables 35b, each of which has the function pattern shown in, for example, Figure 5A. Each of the 
first look-up tables 35b sets the frequencies and the corresponding responses R 1 with respect to the cubic spline inter- 
polating function, which is represented by Formula (39), and for each of a plurality of different image size enlargement 
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scale factors (e.g., enlargement scale factors of 1 .0 t 1 .2, .... 1 .8 Pj and 2.0). The input means 35 is further provided with 
six kinds of second look-up tables 35c, each of which has the function pattern shown in, for example, Figure 5B. Each 
of the second look-up tables 35c sets the frequencies and the corresponding responses R 2 with respect to the B spline 
interpolating function, which is represented by Formula (40), and for each of a plurality of different image size enlarge- 

s ment scale factors (e.g., enlargement scale factors of 1 .0, 1 .2 1 .8, and 2.0). The input means 35 is still further pro- 
vided with a coefficient calculating means 35d. The coefficient calculating means 35d calculates a response of the 
cubic spline interpolating function with respect to an image size enlargement scale factor, which is desired for an inter- 
polation image, and a response R 2 of the B spline interpolating function with respect to the image size enlargement 
scale factor, which is desired for the interpolation image, by making reference to one of the first look-up tables 35b and 

w one of the second look-up tables 35c. which correspond to the image size enlargement scale factor desired for the inter- 
polation image. The coefficient calculating means 35d also calculates the parameter a from an operation carried out 
with Formula (5): 



15 



a=(R-R .,)/(R 2 -R t) (5) 



in accordance with the response R desired for the interpolation image and the calculated responses Ri R 2 . 
The first look-up tables 35b and the second look-up tables 35care stored in a single data base. 
In cases where the input means 35 is constituted in the manner described above, it becomes possible to specify 
the response R, with which the extent of the change in the sharpness can be grasped as an actual feeling. 

20 The image size enlargement scale factor, for which the response of the cubic spline interpolating function and 
the response R 2 of the B spline interpolating function are calculated, is a desired image size enlargement scale factor, 
which is inputted from the input means (not shown) into the multi -formatter 20. 

The information representing the parameter a, which has been inputted directly or has been calculated from the 
inputted response R, is fed into the interpolation coefficient operation means 33. 

25 The interpolation coefficient operation means 33 calculates the new interpolation coefficients a k . 1( a k , a k+1 . and 
a k+2 for each of the values of t according to the value of the parameter a from the received B spline interpolation coef- 
ficients and the received cubic spline interpolation coefficients, which are for each of the values of t, and the parameter 
a. The calculations are carried out with Formulas (41), (42), (43), and (44). 

The information representing the calculated new interpolation coefficients a k . 1t a k , a k+1 , and a k+2 is fed into the 

30 interpolating operation means 34. 

The interpolating operation means 34 calculates the value of the interpolated image signal component Y p , which 
corresponds to the interpolation point Xp, with Formula (38), which is stored therein and serves as the operation formula 
with the third-order spline interpolating function, in accordance with the interpolation coefficients a k ,-, , a k , a k+1 , and a^, 
which have been received from the interpolation coefficient operation means 33, and the original image signal compo- 

35 nents Y k . 1( Y k , Y k+1 , and Y k+2 , which have been received from the image signal storing device 10. 

The thus obtained interpolation image signal S* containing the interpolated image signal components, which corre- 
spond to all interpolation points, is fed into the image reproducing means 40. 

The image reproducing means 40 reproduces a visible image from the received interpolation image signal S\ The 
sharpness of the reproduced visible image can be adjusted easily by merely changing the value of the inputted param- 

40 eter a. In cases where a negative value less than 0 is inputted as the value of the parameter a, an image having a 
sharpness higher than the sharpness of the secondary image, which is obtained with the ordinary cubic spline interpo- 
lating operation, can be obtained. In cases where a value larger than 1 is inputted as the value of the parameter a, a 
smooth image having a sharpness lower than the sharpness of the secondary image, which is obtained with the ordi- 
nary B spline interpolating operation, can be obtained. Also, in cases where the parameter a is set to be equal to a value 

45 larger than 0 and smaller than 1, an image having a sharpness, which is intermediate between the sharpness of the 
secondary image obtained with the cubic spline interpolating operation and the sharpness of the secondary- image 
obtained with the B spline interpolating operation, can be obtained. In cases where the input means 35 has the consti- 
tution for inputting the desired response R, the desired response R can be changed and inputted into the response input 
means 35. In this manner, the value of the parameter a can be changed, and the adjustment of the sharpness can be 

so carried out easily. 

In the image reproducing system provided with this embodiment, the interpolating operation apparatus 30 utilizes 
the primary image signal having been previously stored in the image signal storing device 10. However, the interpolat- 
ing operation apparatus in accordance with the present invention is not limited to this embodiment. For example, an 
image signal representing an image, which has been detected by using an image read-out apparatus illustrated in Fig- 
55 ure 2, may be utilized. 

Specifically, with the image read-out apparatus illustrated in Figure 2, an X-ray image of an object, which has been 
stored on a stimulable phosphor sheet 100, is read from the stimulable phosphor sheet 100. 

The stimulable phosphor sheet 100. on which the X-ray image has been stored, is placed at a predetermined posi- 
tion in a read-out means 50 of an X-ray image read-out apparatus. The stimulable phosphor sheet 100 is then conveyed 
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in a sub-scanning direction, which is indicated by the arrow, Y, by a$ endless belt 52, which is operated by a motor 51 . 
A laser beam 54, which serves as stimulating rays, is produced by a laser beam source 53. The laser beam 54 is 
reflected and deflected by a rotating polygon mirror 56, which is quickly rotated by a motor 55 in the direction indicated 
by the arrow. The laser beam 54 then passes through a converging lens 57, which may be constituted of an f8 lens, or 
the like. The direction of the optical path of the laser beam 54 is then changed by a mirror 58, and the laser beam 54 
impinges upon the stimulable phosphor sheet 100 and scans it in a main scanning direction indicated by the arrow X, 
which direction is approximately normal to the sub-scanning direction indicated by the arrow Y When the stimulable 
phosphor sheet 100 is exposed to the laser beam 54, the exposed portion of the stimulable phosphor sheet 100 emits 
light 59 in an amount proportional to the amount of energy stored thereon during its exposure to the X-rays. The emitted 
light 59 is guided by a light guide member 60 and photoelectrical^ detected by a photomultiplier 61 . 

The light guide member 60 is made from a light guiding material, such as an acrylic plate. The light guide member 
60 has a linear light input face 60a. which is positioned so that it extends along the main scanning line on the stimulable 
phosphor sheet 100. and a ring-shaped light output face 60b, which is positioned so that it is in close contact with a light 
receiving face of the photomultiplier 61. The emitted light 59, which has entered the light guide member 60 at its light 
input face 60a, is guided through repeated total reflection inside of the light guide member 60, emanates from the light 
output face 60b, and is received by the photomultiplier 61. In this manner, the amount of the emitted light 59, which 
amount represents the X-ray image, is converted into an electric signal by the photomultiplier 61 . 

An analog output signal S is generated by the photomultiplier 61 . The analog output signal S is logarithmically 
amplified by a logarithmic amplifier 62, and digitized by an analog-to-digital converter 63. In this manner, an original (pri- 
mary) image signal Dorg is obtained. The original image signal Dorg is fed into the aforesaid multi -formatter 20. 

As described above, the primary image signal, which is utilized in the interpolating operation apparatus 30 in 
accordance with the present invention may be the one, which has been previously stored in the image signal storing 
device 10, or may be the one, which is obtained by reading out an image with the image read-out apparatus illustrated 
in Figure 2. 

Embodiments of spline interpolating operation apparatuses for carrying out the spline interpolating operation meth- 
ods for an image signal in accordance with the present invention will be described hereinbelow. 

Figure 6 is a schematic block diagram showing an image reproducing system provided with an interpolating oper- 
ation apparatus 130, which is an embodiment of the spline interpolating operation apparatus for carrying out the spline 
interpolating operation method for an image signal in accordance with the present invention. With reference to Figure 

6, the image reproducing system comprises an image signal storing device 10 for storing an image signal Sorg repre- 
senting an image, and a multi -formatter 1 20 for carrying out predetermined signal processing on the image signal (here- 
inbelow referred to as the primary image signal or the original image signal) Sorg, which is received from the image 
signal storing device 10, such that an image signal conforming to a predetermined reproduction format may be 
obtained. The image reproducing system also comprises an image reproducing means 40 for reproducing a visible 
image in the predetermined reproduction format from an image signal (hereinbelow referred to as the secondary image 
signal or the interpolation image signal) S\ which has been obtained from the predetermined signal processing carried 
out by the multi-formatter 1 20. 

The multi-formatter 120 carries out the signal processing on the primary image signal Sorg, such that an image sig- 
nal conforming to one of various image reproduction formats may be obtained. Examples of the image reproduction for- 
mats include a format, in which the entire area of a single sheet of film is divided into four different small regions and 
four different images having reduced image sizes are printed respectively in the four regions, a format, in which a single 
large image is printed directly on a single sheet of film, and a format, in which a portion of an image is enlarged and the 
enlarged image portion is printed on a sheet of film. The multi-formatter 120 is provided with the interpolating operation 
apparatus 130 in accordance with the present invention, which carries out an interpolating operation for obtaining the 
secondary image signal (the interpolation image signal) made up of a number of image signal components different 
from that of the image signal components of the primary image signal Sorg when the image size is to be enlarged or 
reduced. 

As illustrated in Figure 7, the primary image signal Sorg used in this embodiment is made up of a series of digital 
image signal components Y k . 2 . Y k . 1( Y k , Y k+1i Y k+2 , ... respectively corresponding to sampling points (picture elements) 
Xk-2- • Xk» *k+i • *k+2 which are sampled with a period of an equal interval and arrayed in one direction. In Figure 

7, the digital image signal components Y M , Y k , Y k+1 , and Y k+2 are indicated by white dots. 

The interpolating operation apparatus 130, which is incorporated in the multi-formatter 120, comprises a storage 
means 131, an input means 132, an interpolation coefficient operation means 133. and an interpolating operation 
means 134. 

The storage means 131 stores information representing interpolation coefficients a^, a k , a^, and a k+2 . which 
respectively correspond to the original image signal components Y k . 1( Y k , Y k+1 , and Y k+2 , in Formula (6). Formula (6) 
serves as a third-order spline interpolating operation formula and represents a secondary image signal component Y p 1 
corresponding to an interpolation point Xp. which is located between original sampling points X k and X k+1 . The interpo- 
lation coefficients a k . 1 , a k , a^, and a»< +2 are stored as being the ones represented by Formulas (13), (14), (15), and 
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(16). 



10 



so 



Y p =a k-1 ^ k-1 +a k Y k +a k+1 Y k+1 +a k+2 Y k+2 ( 6 ) 

a k . 1= -(a/2)t(M) 2 (13) 

a k =(2-a/2)t 3 -(3-a/2)t 2 +1 (14) 

a k+1 =(a/2-2)t 3 +(3- a)t 2 +(a/2)t (15) 

a k+2 =(a/2)t 2 (M) (16) 



in which a represents the parameter determining the sharpness of the secondary image represented by the interpola- 
tion image signal, and t, where 0^t<1 , represents the position of the interpolation point Xp, the position being taken with 
is respect to the picture element X k , that serves as a reference position, and in the direction heading toward the picture 
element X^-j, the lattice interval of the original image signal components being set to be equal to 1 . 

The input means 132 is used to input the arbitrary parameter a, which determines the sharpness of the secondary 
image reproduced from the interpolation image signal. 

The interpolation coefficient operation means 133 calculates the interpolation coefficients a k _ 1( a^ a^-j, and a k+2 
20 in accordance with the parameter a, the calculation being carried out from the interpolation coefficients, which are 
stored in the storage means 131, and the parameter a inputted from the input means 132. 

The interpolating operation means 134 stores Formula (6) shown above as the operation formula with the third- 
order spline interpolating function and calculates the value of the interpolated image signal component Y p , which cor- 
responds to the interpolation point X p( with Formula (6) in accordance with the interpolation coefficients a k . 1( a k , a^, 
25 and ak+2 according to the parameter a, which have been calculated by the interpolation coefficient operation means 
133, and the original image signal components Y k . 1( Y k , Y k+1 , and Y k+2 . 

The interpolation coefficients a^, a^ a k+1 , and a k+2 , which are stored in the storage means 131, are calculated 
previously with the algorithms described below. 

Specifically, as illustrated in Figure 7, the image signal components (the original image signal components), which 

30 represent a series of picture elements X^, X k . 1 , X k , X k+1 , X|< +2 are respectively represented by Y k _ 2 , Y^ , Y k , Y k+1 , 

Y k+2 . ... The third-order spline interpolating function is set for each of sections X^-X^, X k . 1 -X k , X k ~X k+1 , and 
X k+1 -X k+2 . The spline interpolating functions corresponding to the respective sections are represented by f k . 2 , f k _ 1t f k , 
f k+ i, and f k+2 . The interpolating functions are the third-order functions, in which the position in each section serves as 
a variable. 

35 Firstly, in cases where the point taken for interpolation (i.e. the interpolation point) Xp falls within the section 
Xk~Xk+i. the spline interpolating function f k corresponding to the section X k ~X k+1 is represented by Formula (29), 

f k (x)=A k x 3 +B k x 2 +C k x + D k (29) 

40 In the operation formula, Formula (29), the conditions should be satisfied in that the spline interpolating function f k 
should pass through the original sampling points (picture elements), and in that the first-order differential coefficient of 
the spline interpolating function f k should be continuous between adjacent sections. From such conditions, Formulas 
(7), (8), (9), and (10) are derived. 

*k(X k )=Y k (7) 

*k(X k +i)=Y k+1 *(8) 

V(Xk)=W<X k > (9) 
V(X k+1 )=f k+1 '(X k+1 ) (10) 



2 

In these formulas, f k represents the first-order differentiation (3A k x +2B k x+C k ) of the function f k . 
Also, the condition should be satisfied in that the first-order differential coefficient of the spline interpolating function 
55 f k at the picture element X|< should have a predetermined inclination a with respect to the gradient 
(Y k+rYk-iMXk+i'X^) of the original image signal components Y k _i and Y k+1 representing the picture elements X^ 
! and X k+1 , which are located before and after the picture element >4- Therefore, it is necessary for Formula (11) to be 
satisfied: 
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f K '(X K )= a(Y k+1 - Y ^ )/(X ^ -X k ,) (11) 

Further, the condition should be satisfied in that the first-order differential coefficient of the spline interpolating func- 
tion f k at the picture element Xj< +1 should have the predetermined inclination a with respect to the gradient 
(Yk+a'YkWXk+a'Xk) of the ordinal image signal components Y k and Y k+2 representing the picture elements Xk and 
Xk+2. which are located before and after the picture element X^. Therefore, it is necessary for Formula (12) to be sat- 
isfied: 

V(X kt1 )=a(Y k+2 -Y k )/(X k+2 -X k ) (12) 

The interval (i.e., the lattice interval) of each of sections X^-X^, X^-X^ X k ~X k+1 , and Xk +1 ~X k+ 2 is herein set 
to be equal to 1 . Also, the position of the interpolation point Xp. which is taken from the picture element X|< in the direc- 
tion heading toward the picture element X kt1 , is represented by t (0st<1). In such cases, from Formulas (7). (8). (9), 
(10), (1 1). and (12), the formulas shown below obtain. 

f k (0)=D k =Y k 

f k (1)=A k+ B k+ C k +D k =Y k+1 

f k '(0)=C k =a(Y k+1 -Y k . 1 )/2 

f k '(1 )=3A k +2B k+ C k =a(Y k+2 -Y k )/2 

Therefore, the formulas shown below obtain. 

A k =(a/2)Y k+2+ (a/2-2)Y k+1 +(2-a/2)Y k -(a/2)Y k ., 

B k =-(a/2)Y k+2 +(3-a)Y k+1 +(a/2-3)Y k+ aY M 

C k =(a/2)Y k+l -(a/2)Y k . 1 

D k =Y k 

As described above, the variable conversion of X=t is carried out, and therefore the spline interpolating function 
f k (x) is represented by the formula shown below. 

f K M=f k (t) 

Therefore, an interpolated image signal component Y p corresponding to the interpolation point X p may be represented 
by Formula (30). 

Y P=Mt)=A K t 3 +B k t 2 + C k t+D k (30) 
Substituting the coefficients \, B|<. C*. and D k into Formula (30) yields 

Y p ={(a/2)Y k+2+ (o/2-2)Y k+1 +(2-a/2)Y k -(a/2)Y k .,}t 3 +{-(a/2)Y k+2 +(3-a)Y k+1 

+(a/2-3)Y k +aY }t 2 +{(a/2)Y k+1 -(a/2)Y k . 1 }t+Y k 

Arranging this formula with respect to the image signal components Y k .i. Y k . Y k+1 , and Y k+2 yields Formula (45). 

Y_ = {at (M) 2 12) Y M + {(2-o/2) t 3 + (a/2-3) t 2 +1} Y k 

(45) 

+ {(a/2-2) t 3 + (3-a) t 2 +at/2} Y k+1 + {-at 2 (t-1) 12) Y k+2 

The interpolation coefficients a^. a^, ai< + i, and a k+2 corresponding to the original image signal components Y^. 
Y k+1 , and Y k+2 are calculated with the algorithms described above. 
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The interpolation coefficient operation means 133 calculates the values of the interpolation coefficients a^. a k , 
ak +1 , and a k+2 . which coefficients have been stored in the storage means 131, in accordance with an arbitrary param- 
eter a, which corresponds to a desired level of sharpness and has been inputted from the input means 132. 

As will be understood from Formulas (1 1 ) and (12) and Figure 7, as the value of the parameter a is set to be larger, 
5 the contrast of the secondary image becomes higher, and an image having a higher sharpness is obtained. As the 
value of the parameter a is set to be smaller, the contrast of the secondary image becomes lower, and a smoother 
image having a lower sharpness is obtained. 

An actual image is composed of the picture elements arrayed in two-dimensional directions. Therefore, the param- 
eter t, which is used in the aforesaid interpolation coefficient formulas, may be replaced by t x with respect to one of the 
w two directions in the array and may be replaced by t y with respect to the other direction in the array. The interpolation 
coefficients may thus be determined for each of the directions in the array. In this specification, as an aid in facilitating 
the explanation, the operation is described with respect to only one of the two directions in the array of the picture ele- 
ments. 

How the image reproducing system provided with this embodiment operates will be described hereinbelow. 
15 Firstly, the multi-formatter 120 reads the primary image signal Sorg from the image signal storing device 10. Also, 
in order to obtain the secondary image signal, which represents an enlarged image corresponding to the image size 
enlargement scale factor inputted from the input means (not shown), the multi-formatter 120 feeds the primary image 
signal Sorg into the interpolating operation apparatus 130 incorporated in the multi-formatter 120. 

In the interpolating operation apparatus 130, the primary image signal Sorg is fed into the interpolating operation 
20 means 134. 

The storage means 131 sets the values of t in Formulas (13), (14), (15), and (16) in accordance with the image size 
enlargement scale factor having been inputted from a different input means (not shown). For example, in cases where 
an image size enlargement scale factor of 2 is inputted, values of 0.0 and 0.5 are set as the values of t. In cases where 
an image size enlargement scale factor of 4 is inputted, values of 0.0, 0.25, 0.5, and 0.75 are set as the values of t. In 
25 cases where an image size enlargement scale factor of 10 is inputted, values of 0.0, 0.1 , 0.2 0.9 are set as the val- 
ues of t. Information representing the interpolation coefficients a k . t , a k , a K+1 , and a k+2 , which are represented by For- 
mulas (13), (14), (15), and (16) and are for each of the thus set values of t, is fed into the interpolation coefficient 
operation means 133. 

Also, information representing a value of the parameter a corresponding to a level of sharpness desired for the sec- 
30 ondary image is inputted into the input means 132. The information representing the value of the parameter a is also 
fed into the interpolation coefficient operation means 133. 

The interpolation coefficient operation means 133 calculates the interpolation coefficients a k . 1( a^ a k+1 , and a k+2 
for each of the values of t according to the value of the parameter a from the received interpolation coefficients a k _j . a k . 
ak+1. and ak+2. which are for each of the values of t, and the parameter a. The calculations are carried out with Formulas 
35 (13). (14), (15), and (16). 

The information representing the calculated interpolation coefficients a^, a^ a^, and a k+2 is fed into the interpo- 
lating operation means 134. 

The interpolating operation means 134 calculates the value of the interpolated image signal component Yp, which 
corresponds to the interpolation point Xp, with Formula (6), which is stored therein and serves as the operation formula 

40 with the third-order spline interpolating function, in accordance with the interpolation coefficients a^ , a k , a k+1 , and a k+2 , 
which have been received from the interpolation coefficient operation means 133, and the original image signal com- 
ponents Y k _ 1f Y k , Y k+1 , and Y k+2 , which have been received from the image signal storing device 10. 

The thus obtained interpolation image signal S* containing the interpolated image signal components, which corre- 
spond to all interpolation points, is fed into the image reproducing means 40. 

45 The image reproducing means 40 reproduces a visible image from the received interpolation image signal S*. The 
sharpness of the reproduced visible image can be adjusted easily by merely changing the value of the inputted param- 
eter a. Specifically, with this embodiment of the spline interpolating operation apparatus for an image signal in accord- 
ance with the present invention, wherein two or more kinds of spline interpolating operations are not combined with one 
another, the sharpness can be adjusted easily. Also, the value of the parameter a can be set freely, and therefore the 

so degree of freedom of the sharpness adjustment can be kept high. 

Figures 8A, 8B, and 8C show the graphs of the interpolation coefficients obtained when the value of the parameter 
a is set to be equal to 0.4, 1 .0, and 1 .6. respectively. 

In the image reproducing system provided with this embodiment of the spline interpolating operation apparatus for 
an image signal in accordance with the present invention, the interpolating operation apparatus 130 utilizes the primary 

55 image signal having been previously stored in the image signal storing device 10. However, the spline interpolating 
operation apparatus in accordance with the present invention is not limited to this embodiment. For example, an image 
signal representing an image, which has been detected by using an image read-out apparatus illustrated in Figure 2, 
may be utilized. 

As described above, the primary image signal, which is utilized in the interpolating operation apparatus 130 in 
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accordance with the present invention may be the one, which has peen previously stored in the image signal storing 
device 10, or may be the one, which is obtained by reading out an image with the image read-out apparatus illustrated 
in Figure 2. 

Also, in this embodiment of the spline interpolating operation apparatus for an image signal in accordance with the 
present invention, the information representing the interpolation coefficients, which are represented by Formulas (13), 
(14), (15), and (16), are stored in the storage means 131. Alternatively, as the interpolation coefficients, for example, 
those represented by Formulas (21), (22), (23), and (24) shown below may be employed. 



a k . 1= -(p+1/2)t 3 +(3p/2+1)t 2 -(1/2)t-p/2 (21) 

a k =3(p+1/2)t 3 -(9p/2+5/2)t 2 +P+1 (22) 

a k+1 =-3(p+1/2)t 3 +(9p/2+2)t 2 +(1/2)t-p/2 (23) 

a k+2 =(p+1/2)t 3 -(3p/2+1/2)t 2 (24) 



In Formulas (21), (22), (23), and (24), p represents an arbitrary parameter, which determines the sharpness of the 
secondary image represented by the interpolation image signal. In lieu of the aforesaid parameter a, the parameter p 
is inputted from the input means 132. The interpolation coefficients a^, a^ a^, and a^ are calculated previously 
with the algorithms described below. 

Specifically, as illustrated in Figure 7, the image signal components (the original image signal components), which 

represent a series of picture elements X^. X^.,, X^ Xk +1 , are respectively represented by Y k _ 2 , Y k . 1( Y k , Y k+1 , 

Y k+2 , ... The third-order spline interpolating function is set for each of sections X k . 2 -X k . 1 , X^-X^ X^-X^-,, and 
x k+i~ x k+2- The spline interpolating functions corresponding to the respective sections are represented by f k _ 2 , f k -i. f k . 
f k+1 , and f k+2 . The interpolating functions are the third-order functions, in which the position in each section serves as 
a variable. 

Firstly, in cases where the point taken for interpolation (i.e. the interpolation point) Xp falls within the section 
Xk~ x k+1. the spline interpolating function f k corresponding to the section X k -X k+1 is represented by Formula (29). 

f k (x)=A k x 3 +B k x 2 +C k x+D k (29) 

In the operation formula, Formula (29), the conditions should be satisfied in that the spline interpolating function f k 
should slightly deviate from the original sampling points (picture elements) in order to carry out an smooth interpolating 
operation, and in that a non-uniform density portion should not occur in the interpolation image, which is obtained from 
the interpolating operation carried out on an original image having uniform density distribution. From such conditions, 
Formulas (1 7) and (18) are derived. 

f k (X k )=-0.5pY ^ +(1+p)Y k -0.5pY k+1 (17) 
* k(* k + i )=-0-5p Y k+ (1 +p)Y k+1 -0.5PY k+2 (1 8) 

The parameter p in Formulas (17) and (18) are the parameter, which determines the sharpness of the secondary 
image represented by the interpolation image signal. Formulas (17) and (18) mean that the secondary image signal 
components corresponding to the points X k and at both ends of the section X k -X k+1 are raised or lowered by the 
inputted value of p from the original image signal components, and that halves of p are allocated to the adjacent points 
on both sides and each of the original image signal components on both sides is subtracted by the amount of contribu- 
tion of the allocation. 

From the condition in that the first-order differential coefficient of the spline interpolating function f k should be con- 
tinuous between adjacent sections, Formulas (9) and (10) are derived. 

w k )=f k y(x k ) o) 

V<x k+l )=W(x k+1 ) do) 

Also, the condition should be satisfied in that the first-order differential coefficient of the spline interpolating function 
f k at the picture element ^ should coincide with the gradient (Y k+1 -Y ko )/(X k+1 -X ko ) of the original image signal com- 
ponents Y k-1 and Y k+1 representing the picture elements X^ and X k+1 , which are located before and after the picture 
element X^ Therefore, it is necessary for Formula (19) to be satisfied. 
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^'(XhWY^-Y^WX^-X^) (19) 

Further, the condition should be satisfied in that the first-order differential coefficient of the spline interpolating func- 
tion f k at the picture element should coincide with the gradient (Y k+2 -Y k )/(X k+2 -X k ) of the original image signal 
components Y k and Y k+2 representing the picture elements Xj< and X k+2 , which are located before and after the picture 
element X k+1 . Therefore, it is necessary for Formula (20) to be satisfied. 

*K*(X k+1 )=(Y k+2 -Y K )/(X k+2 -X k ) (20) 

The interval (i.e., the lattice interval) of each of sections X^-X^, X^-X^ X k -X k+1 , and X k+1 ~X k+2 is herein set 
to be equal to 1 . Also, the position of the interpolation point Xp, which is taken from the picture element X k in the direc- 
tion heading toward the picture element X k+1( is represented by t (0^t<1). In such cases, from Formulas (9), (10), (17), 
(18), (19), and (20), the formulas shown below obtain. 

f k (0)=D k =-0.5pY k . 1+ (1+p)Y k -0.5pY k+l 

f k (1)=A k+ B k+ C k+ D k 
=-0.5pY k +(1+p)Y k+1 -0.5pY 

k+2 

f k '(0)«C k -(Y k+1 -Y t1 )y2 
V(1)=3A k +2B k +C k =(Y k+2 -Y k )/2 
Therefore, the formulas shown below obtain. 

A k =(1/2+p)Y k+2 -(3/2+3p)Y k+1 +(3p+3/2)Y k -(P+1/2)Y k ., 
B k =-(1/2+3p/2)Y k+2 +(9p/2+2)Y k+1 -(9p/2+5/2)Y k +(1+3p/2)Y k . 1 

C k =(1/2)Y k+1 -(1/2)Y k .., 
DfcHP/SJY^+O+pjY^p^Y^ 

As described above, the variable conversion of X=t is carried out, and therefore the spline interpolating function 
f k (x) is represented by the formula shown below. 

Mx)=f k (t) 

Therefore, an interpolated image signal component Y p corresponding to the interpolation point Xp may be represented 
by Formula (30). 

Y p =f k (t)=A k t 3 + B k t 2 +C k t+D 1< (30) 
Substituting the coefficients A k , B k , C^, and D k into Formula (30) yields 

Y p ={(1/2+p)Y k+2 -(3/2+3p)Y k+l +(3p+3/2)Y k -(p+1/2)Y k . 1 }t 3 +{-(1/2+3 p/2)Y k+2 +(9p/2+2)Y k+1 
-OP^+^Y^I+SP^., }t 2 +{(1/2)Y k+1 -(1/2)Y k .., }t+KP/2)Y k+1 +(1+p)Y k -(p/2)Y 

Arranging this formula with respect to the image signal components Y^. Y k , Y k+1 , and Y k+2 yields Formula (46). 
Y„ = {-(P+1/2) t 3 + (3P/2+1) t 2 - (1/2) t-p/2} Y k ., + {3 (P+1/2) t 3 - (9p/2+5/2) t 2 + p+1} Y k 

(46) 

+ {-3 (P+1/2) t 3 + (9P/2+2) t 2 + (1/2) t-p/2} Y k+1 + {(P+1/2) t 3 - (3p/2+1/2)t 2 )Y k+2 
The interpolation coefficients , ak. ak+i , and corresponding to the original image signal components , 
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Y k , Y k+1 , and Y k+2 are calculated with the algorithms described abovp. 

In such cases, the interpolation coefficient operation means 133 calculates the interpolation coefficients a^, a k , 
a k+1 , and a k+2 in accordance with the parameter p, the calculation being carried out from the interpolation coefficients, 
which are stored in the storage means 131 , and the parameter p inputted from the input means 132. The interpolating 
operation means 134 calculates the value of the interpolated image signal component Y p , which corresponds to the 
interpolation point X p , with Formula (6) in accordance with the interpolation coefficients a k _ 1( a k , a k+1 , and a k+2 accord- 
ing to the parameter p, which have been calculated by the interpolation coefficient operation means 133, and the orig- 
inal image signal components Y k . 1( Y k , Y k+1 , and Y k+2 . 

Figures 9A, 9B, and 9C show the graphs of the interpolation coefficients obtained when the value of the parameter 
p is set to be equal to -0.5. 0, and 0.5, respectively: As will be understood from Figures 9A, 9B, and 9C, as the value of 
the parameter p is set to be larger, the contrast of the secondary image becomes higher, and an image having a higher 
sharpness is obtained. As the value of the parameter p is set to be smaller, the contrast of the secondary image 
becomes lower, and a smoother image having a lower sharpness is obtained. 

As the interpolation coefficients a k . 1 , a k , a k+1 , and a k+2 stored in the storage means 131, those obtained from com- 
binations of the aforesaid parameters a and p may be employed. Specifically, the mean value of the values of Formulas 
(13) and (21), which is calculated with Formula (25), may be employed as the interpolation coefficient a k .-|. The mean 
value of the values of Formulas (14) and (22), which is calculated with Formula (26), may be employed as the interpo- 
lation coefficient a k . The mean value of the values of Formulas (15) and (23), which is calculated with Formula (27), may 
be employed as the interpolation coefficient a k+1 . Also, the mean value of the values of Formulas (16) and (24), which 
is calculated with Formula (28), may be employed as the interpolation coefficient a k+2 . 

a k . 1 =-{(a+2p+1)/4}t 3 +{(2a+3p+2)/4}t 2 -{(a+1)/4}t-p/4 (25) 

a k ={(-a+6p+7)/4}t 3 +{(a-9p-11)/4}t 2 +p/2+1 (26) 

a k+1 ={(a-6p-7)/4}t 3 +{(-2a+9p+10)/4}t 2 +{(a+1 )/4}t-p/4 (27) 

a k+2 -{(a + 2p + 1)/4}t 3 -{(a+3p+1)/4}t 2 (28) 

In such cases, the interpolation coefficient operation means 133 calculates the interpolation coefficients a k . 1 , a k , 
a^, and a|< +2 in accordance with the parameters a and p, the calculation being carried out from the interpolation coef- 
ficients, which are stored in the storage means 131, and the parameters a and p inputted from the input means 132. 
The interpolating operation means 134 calculates the value of the interpolated image signal component Y p , which cor- 
responds to the interpolation point Xp, with Formula (6) in accordance with the interpolation coefficients a^-i, a k( a^, 
and a k+2 according to the parameters a and p, which coefficients have been calculated by the interpolation coefficient 
operation means 133, and the original image signal components Y k . 1t Y k , Y k+1 , and Y k+2 . 

Figures 10A, 10B t and 10C are graphs showing interpolation coefficients obtained when the values of the param- 
eters (a, p) are set to be equal to (0.4, -0.5), (1 .0, 0.0), and (1 .6, 0.5), respectively. As will be understood from Figures 
10A, 10B, and 10C, as the value of the parameter a is set to be larger, or as the value of the parameter p is set to be 
larger, the contrast of the secondary image becomes higher, and an image having a higher sharpness is obtained. As 
the value of the parameter a is set to be smaller, or as the value of the parameter p is set to be smaller, the contrast of 
the secondary image becomes lower, and a smoother image having a lower sharpness is obtained. 

Claims 

1. An interpolating operation method for an image signal, comprising the steps of: 

i) obtaining an original image signal, which represents an original image and is made up of a series of original 
image signal components Yij, 

ii) linearly combining interpolation coefficients Bi] and Cij, which correspond to each other and are set for each 
of the original image signal components Yij, in two different interpolating functions f and g for obtaining two 
interpolation images having different levels of sharpness, which functions are represented by Formulas (1 ) and 
(2), said linear combination being carried out with Formula (3), a new interpolation coefficient Aij being 
obtained from said linear combination, and 

iii) carrying out an interpolating operation on the original image signal components Yij by using an interpolating 
function h having the new interpolation coefficient Aij, which function is represented by Formula (4), an inter- 
polation image signal being obtained from said interpolating operation, said interpolation image signal being 
made up of a series of image signal components, which occur at intervals different from those of the original 
image signal components Yij. 
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f=IBij.Yij a (1) 

g=£Cij-Yij (2) 
5 Aij=(1-a)Bij+aCij (3) 

h=LAij • Yij (4) 

in which i=1 , 2 and j=1 , 2 

w wherein the coefficient a in Formula (3) is set to be one of real numbers including a range smaller than 

0 and/or a range larger than 1 . 

2. An interpolating operation method for an image signal as defined in Claim 1 wherein: 

is a) a plurality of first look-up tables are prepared, each of said first look-up tables setting the frequencies and 

the corresponding responses with respect to one of said two interpolating functions and for each of a plu- 
rality of different image size enlargement scale factors, 

b) a plurality of second look-up tables are prepared, each of said second look-up tables setting the frequencies 
and the corresponding responses R 2 with respect to the other interpolating function and for each of a plurality 

20 of different image size enlargement scale factors, 

c) a response R 1 of said one interpolating function with respect to an image size enlargement scale factor, 
which is desired for an interpolation image, and a response R 2 of said other interpolating function with respect 
to said image size enlargement scale factor, which is desired for the interpolation image, are calculated by 
making reference to a first look-up table and a second look-up table, which correspond to said image size 

25 enlargement scale factor desired for the interpolation image, and 

d) said coefficient a is calculated from an operation carried out with Formula (5): 



ct=(R-R .,)/(R 2 -R <,) (5) 

30 in accordance with a response R desired for the interpolation image, the calculated response R-i of said one 

interpolating function, and the calculated response R 2 of said other interpolating function. 

3. An interpolating operation method for an image signal as defined in Claim 1 wherein one of the two different inter- 
polating functions for obtaining two interpolation images having different levels of sharpness is a B spline interpo- 
35 lating operation function, and the other interpolating function is a cubic spline interpolating operation function. 



4. An interpolating operation method for an image signal as defined in Claim 2 wherein one of the two different inter- 
polating functions for obtaining two interpolation images having different levels of sharpness is a B spline interpo- 
lating operation function, and the other interpolating function is a cubic spline interpolating operation function. 

40 

5. An interpolating operation method for an image signal as defined in Claim 1 wherein the original image is a radia- 
tion image. 

6. An interpolating operation method for an image signal as defined in Claim 5 wherein the radiation image has been 
45 stored on a stimulable phosphor sheet. 

7. An interpolating operation method for an image signal as defined in Claim 6 wherein the original image signal rep- 
resenting the radiation image is obtained by exposing the stimulable phosphor sheet to stimulating rays, which 
cause the stimulable phosphor sheet to emit light in proportion to the amount of energy stored thereon during its 

so exposure to the radiation, and photoelectrical^ detecting the emitted light. 

8. An interpolating operation method for an image signal as defined in Claim 1 wherein the original image has been 
recorded on photographic film. 

55 9. An interpolating operation apparatus for an image signal, wherein: 

i) an original image signal, which represents an original image and is made up of a series of original image sig- 
nal components Yij, is obtained, 

ii) interpolation coefficients Bij and Cij, which correspond to each other and are set for each of the original 
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image signal components Yij, in two different interpolating functions f and g for obtaining two interpolation 
images having different levels of sharpness, which functions are represented by Formulas (1) and (2), are lin- 
early combined with each other, said linear combination being carried out with Formula (3), a new interpolation 
coefficient Aij being obtained from said linear combination, and 

iii) an interpolating operation is carried out on the original image signal components Yij by using an interpolat- 
ing function h having the new interpolation coefficient Aij, which function is represented by Formula (4), an 
interpolation image signal being obtained from said interpolating operation, said interpolation image signal 
being made up of a series of image signal components, which occur at intervals different from those of the orig- 
inal image signal components Yij, 



f=£Bij- Yij (1) 

g=£Cij • Yij (2) 

Aij=(1-a)Bij+aCij (3) 

h=£Aij * Yij (4) 



in which i=1, 2 and j=1, 2 

the apparatus comprising: 



1) a storage means for storing information representing the interpolation coefficients Bij and Cij. 

2) an input means for inputting the coefficient a as one of real numbers including a range smaller than 0 
and/or a range larger than 1, the coefficient a determining the sharpness of a secondary image, which is 
reproduced from the interpolation image signal, 

3) an interpolation coefficient operation means for calculating the interpolation coefficient Aij in accord- 
ance with the coefficient a, the calculation being carried out from the interpolation coefficients Bij and Cij, 
which are stored in said storage means, and the coefficient a inputted from said input means, and 

4) an interpolating operation means for storing Formula (4) as the operation formula, and calculating the 
value of an interpolated image signal component Y p , which corresponds to an interpolation point Xp, with 
Formula (4) in accordance with the interpolation coefficient Aij, which has been calculated by said interpo- 
lation coefficient operation means, and each of the original image signal components Yij. 

0. An interpolating operation apparatus for an image signal as defined in Claim 9 wherein said input means is pro- 
vided with: 



a) a response input means for inputting a response R desired for the interpolation image, 

b) a plurality of first look-up tables, each of said first look-up tables setting the frequencies and the correspond- 
ing responses with respect to one of said two interpolating functions and for each of a plurality of different 
image size enlargement scale factors, 

c) a plurality of second look-up tables, each of said second look-up tables setting the frequencies and the cor- 
responding responses R 2 with respect to the other interpolating function and for each of a plurality of different 
image size enlargement scale factors, and 

d) a coeff icient calculating means for calculating a response R 1 of said one interpolating function with respect 
to an image size enlargement scale factor, which is desired for an interpolation image, and a response R 2 of 
said other interpolating function with respect to said image size enlargement scale factor, which is desired for 
the interpolation image, by making reference to a first look-up table and a second look-up table, which corre- 
spond to said image size enlargement scale factor desired for the interpolation image, and 

for calculating said coefficient a from an operation carried out with Formula (5): 

a=(R-R ^/(Rg-R ^ (5) 



in accordance with the response R desired for the interpolation image, the calculated response R 1 of said one 
interpolating function, and the calculated response R 2 of said other interpolating function. 

1. An interpolating operation apparatus for an image signal as defined in Claim 9 wherein one of the two different 
interpolating functions for obtaining two interpolation images having different levels of sharpness is a B spline inter- 
polating operation function, and the other interpolating function is a cubic spline interpolating operation function. 
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12. An interpolating operation apparatus for an image signal as defined in Claim 10 wherein one of the two different 
interpolating functions for obtaining two interpolation images having different levels of sharpness is a B spline inter- 
polating operation function, and the other interpolating function is a cubic spline interpolating operation function. 



5 13. An interpolating operation apparatus for an image signal as defined in Claim 9 wherein the original image is a radi- 
ation image. 

14. An interpolating operation apparatus for an image signal as defined in Claim 13 wherein the radiation image has 
been stored on a stimulable phosphor sheet. 

10 

15. An interpolating operation apparatus for an image signal as defined in Claim 14 wherein the original image signal 
representing the radiation image is obtained by exposing the stimulable phosphor sheet to stimulating rays, which 
cause the stimulable phosphor sheet to emit light in proportion to the amount of energy stored thereon during its 
exposure to the radiation, and photoelectrical^ detecting the emitted light. 

15 

1 6. An interpolating operation apparatus for an image signal as defined in Claim 9 wherein the original image has been 
recorded on photographic film. 



17. A spline interpolating operation method for an image signal, wherein an interpolating operation using a third-order 
20 spline interpolating function is carried out on an original image signal, which represents an original image and is 
made up of a series of original image signal components, and an interpolation image signal is obtained with For- 
mula (6), said interpolation image signal being made up of a series of image signal components, which occur at 
intervals different from those of the original image signal components, 



25 Y p =a k- 1 Y k-1 +a k Y k +a k+ 1 Y k+ 1 +a k+2 Y k+2 ( 6 ) 



in which Y k _ 1t Y k , Y k+1 , and Y k+2 represent the original image signal components representing picture elements X^ 
1, X k , Xk+L and X k+2 in the original image, Y p represents the interpolated image signal component corresponding 
to an interpolation point Xp located between the picture elements X k and X k+1 , and a k . A , a k , a k+1 , and a k+2 repre- 
30 sent the interpolation coefficients, 

the method comprising the steps of: 

i) calculating the interpolation coefficients a^, a k . a k+1 , and a k+2 respectively corresponding to the original 
image signal components Y k-1 , Y k , Y k+1 , and Y k+2 such that, in cases where the original image signal compo- 

35 nents representing the two picture elements X^ and Xj< located before the interpolation point Xp, which is 

located between the picture elements X k and X k+1 in the original image, are represented by Y^ and Y k , the 
original image signal components representing the two picture elements X k+1 and X k+2 located after the inter- 
polation point Xp are represented by Y k+1 and Y k+2 , the third-order spline interpolating function between the 
picture elements and X k+1 is represented by f k , the third-order spline interpolating function between the pic- 

40 ture elements X k .-| and X« is represented by f k . 1( and the third-order spline interpolating function between the 

picture elements Xk +1 and X k+2 is represented by f k+1 : 



a) the spline interpolating function f k at the picture elements X k and X k+1 may satisfy the original image sig- 
nal components Y k and Y k+1 , as represented by Formulas (7) and (8): 

MX k )=Y k (7) 

fk(X k+ i)=V k+ i < 8 > 

so b) the first-order differential coefficient of the spline interpolating function f k at the picture element Xk may 

coincide with the first-order differential coefficient of the spline interpolating function f k .-| at the picture ele- 
ment X|<, as represented by Formula (9): 



V(Xk)=fk.i'<Xk) < 9 > 

c) the first-order differential coefficient of the spline interpolating function f k at the picture element X|< +1 may 
coincide with the first-order differential coefficient of the spline interpolating function f k+1 at the picture ele- 
ment Xk +1t as represented by Formula (10): 



33 



EP 0 753 828 A2 



* k k+lM k+1 k+l)' 



(10) 



d) the first-order differential coefficient of the spline interpolating function f k at the picture element X k may 
have an inclination of an arbitrary parameter a with respect to the gradient of the original image signal 
components Y k _ t and Y k+1 representing the picture elements and X k+1 , which are located before and 
after the picture element X*. as represented by Formula (1 1), the arbitrary parameter a being selected pre- 
viously and determining the sharpness of a secondary image represented by the interpolation image sig- 
nal, 

V(* k)=«(Y k+ i -Y k _! )/(X h+1 -X k0 ) (11) 

and 

e) the first-order differential coefficient of the spline interpolating function f k at the picture element X k+1 
may have an inclination of said parameter a with respect to the gradient of the original image signal com- 
ponents Y k and Y k+2 representing the picture elements X k and X k+2 , which are located before and after the 
picture element X^, as represented by Formula (12): 

f k '(X k+1 )=a(Y k+2 -Y k )/(X k+2 -X k ) (12) 

and 

ii) calculating the interpolated image signal component Y p corresponding to the interpolation point X p in 
accordance with the calculated interpolation coefficients a k . 1f a k , a^, and a^. and the original image signal 
components Y k . 1t Y k , Y k+1 , and Y k+2 . 

18. A spline interpolating operation method for an image signal as defined in Claim 17 wherein the calculated interpo- 
lation coefficients a k . 1p a k+1 , and a k+2 are respectively represented by Formulas (13). (14), (15), and (16): 



a k-1 =-(ct/2)t(t-1) 2 (13) 

a k =(2-a/2)t 3 -(3-a/2)t 2 +1 (14) 

a k+1 =(o/2-2)t 3 +(3-a)t 2 +(a/2)t (15) 

a k+2 =(a/2)t 2 (M) (16) 



in which t, where 0^t<1 , represents the position of the interpolation point Xp, the position being taken with respect 
to the picture element X,< ( that serves as a reference position, and in the direction heading toward the picture ele- 
ment X k+1f the lattice interval of the original image signal components being set to be equal to 1 . 

19. A spline interpolating operation method for an image signal as defined in Claim 17 wherein the original image is a 
radiation image. 

20. A spline interpolating operation method for an image signal as defined in Claim 19 wherein the radiation image has 
been stored on a stimulable phosphor sheet. 

21. A spline interpolating operation method for an image signal as defined in Claim 20 wherein the original image sig- 
nal representing the radiation image is obtained by exposing the stimulable phosphor sheet to stimulating rays, 
which cause the stimulable phosphor sheet to emit light in proportion to the amount of energy stored thereon during 
its exposure to the radiation, and photoelectrical^ detecting the emitted light. 

22. A spline interpolating operation method for an image signal as defined in Claim 17 wherein the original image has 
been recorded on photographic film. 

23. A spline interpolating operation method for an image signal, wherein an interpolating operation using a third-order 
spline interpolating function is carried out on an original image signal, which represents an original image and is 
made up of a series of original image signal components, and an interpolation image signal is obtained with For- 
mula (6), said interpolation image signal being made up of a series of image signal components, which occur at 
intervals different from those of the original image signal components, 
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^ p =a KO Y k ^ +a k Y k +a k+1 Y k + ^ +a k+2 Y k+2 (6) 

in which Y k . 1( Y k , Y k+1 , and Y k+2 represent the original image signal components representing picture elements 
1t Xk, X k+1 , and X k+2 in the original image, Y p represents the interpolated image signal component corresponding 
to an interpolation point X p located between the picture elements Xj< and X k+1 , and a k . 1f a k , a|< +1f and a k+2 repre- 
sent the interpolation coefficients, 

the method comprising the steps of: 

i) calculating the interpolation coefficients a^-,, a k . a k+1 , and a k+2 respectively corresponding to the original 
image signal components Y k . 1t Y k , Y k+1 , and Y k+2 such that, in cases where the original image signal compo- 
nents representing the two picture elements X^ and Xj< located before the interpolation point Xp, which is 
located between the picture elements X k and X|<+1 in the original image, are represented by Y^ and Y k , the 
original image signal components representing the two picture elements X k+1 and Xj< +2 located after the inter- 
polation point Xp are represented by Y k+1 and Y k+2 , the third-order spline interpolating function between the 
picture elements X k and X k+1 is represented by f k , the third-order spline interpolating function between the pic- 
ture elements X k _-, and X k is represented by f k . 1f and the third-order spline interpolating function between the 
picture elements Xj< +1 and X k+2 is represented by f k+1 : 

a) the spline interpolating function f k at the picture elements X k and X k+1 may slightly deviate from the orig- 
inal image signal components Y k and Y k+1 in accordance with an arbitrary parameter p, as represented by 
Formulas (17) and (18), the arbitrary parameter p being selected previously and determining the sharp- 
ness of a secondary image represented by the interpolation image signal, 

f k (X k )=-0.5p Y k ., +(1 +p)Y k -0.5pY k+1 (1 7) 

f k (X k+1 )=-0.5p Y k +(1 +P)Y k+1 -0.5PY k+2 (1 8) 

b) the first-order differential coefficient of the spline interpolating function f k at the picture element X k may 
coincide with the first-order differential coefficient of the spline interpolating function f k .-, at the picture ele- 
ment X k , as represented by Formula (9): 

WkM^W (9) 

c) the first-order differential coefficient of the spline interpolating function f k at the picture element X|< +1 may 
coincide with the first-order differential coefficient of the spline interpolating function f k+1 at the picture ele- 
ment X k+1 , as represented by Formula (10): 

K+1 M k+1 

XX k+1 ) do) 

d) the first-order differential coefficient of the spline interpolating function f k at the picture element Xk may 
coincide with the gradient of the original image signal components Y k-1 and Y k+1 representing the picture 
elements X^ and Xj< +1 , which are located before and after the picture element X^ as represented by For- 
mula (19): 

fkXX.MY^-Y^^X^.-X^) (19) 

and 

e) the first-order differential coefficient of the spline interpolating function f k at the picture element Xj< +1 
may coincide with the gradient of the original image signal components Y k and Y k+2 representing the pic- 
ture elements ^ and X k+2 , which are located before and after the picture element X k+1 , as represented by 
Formula (20) : 

* k*(X k+ i )=(Y k+2 -Y k )/(X k+2 -X k ) (20) 

and 

ii) calculating the interpolated image signal component Y p corresponding to the interpolation point Xp in 
accordance with the calculated interpolation coefficients a^, a*. a k+ i , and a^, and the original image signal 
components Y k . 1( Y k , Y k+1 , and Y k+2 . 
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24. A spline interpolating operation method for an image signal as defined in Claim 23 wherein the calculated interpo- 
lation coefficients a k _ 1( a^ a^, and a k+2 are respectively represented by Formulas (21), (22), (23), and (24): 

a k .,='(p+l/2)\ 3 +(3p/2+1)\ 2 -(1/2)t-p/2 (21) 

a k =3(p+1/2)t 3 -(9p/2+5/2)t 2 +p+1 (22) 

a k+1 =-3(p+1/2)t 3 +(9p/2+2)t 2 +(1/2)t-p/2 (23) 

a k+2 =(P+1/2)t 3 -(3p/2+1/2)t 2 (24) 

in which t, where 0^t<1 , represents the position of the interpolation point Xp, the position being taken with respect 
to the picture element X k , that serves as a reference position, and in the direction heading toward the picture ele- 
ment Xj^-i, the lattice interval of the original image signal components being set to be equal to 1. 

25. A spline interpolating operation method for an image signal as defined in Claim 23 wherein the original image is a 
radiation image. 

26. A spline interpolating operation method for an image signal as defined in Claim 25 wherein the radiation image has 
been stored on a stimulable phosphor sheet. 

27. A spline interpolating operation method for an image signal as defined in Claim 26 wherein the original image sig- 
nal representing the radiation image is obtained by exposing the stimulable phosphor sheet to stimulating rays, 
which cause the stimulable phosphor sheet to emit light in proportion to the amount of energy stored thereon during 
its exposure to the radiation, and photoelectrical^ detecting the emitted light. 

28. A spline interpolating operation method for an image signal as defined in Claim 23 wherein the original image has 
been recorded on photographic film. 

29. A spline interpolating operation method for an image signal, wherein an interpolating operation using a third-order 
spline interpolating function is carried out on an original image signal, which represents an original image and is 
made up of a series of original image signal components, and an interpolation image signal is obtained with For- 
mula (6), said interpolation image signal being made up of a series of image signal components, which occur at 
intervals different from those of the original image signal components, 

Y P = a k- 1 Y k-1 +a k Y k+ a k + i Y k+1 +a k+2 Y k+2 (6) 

in which Y k . 1( Y k , Y k+1 , and Y k+2 represent the original image signal components representing picture elements 
i. Xk. X k+1. and *k+2 in the original image, Y p represents the interpolated image signal component corresponding 
to an interpolation point Xp located between the picture elements Xk and Xk +1> and a^, a k( a k+1 , and a k+2 repre- 
sent the interpolation coefficients, 

the method comprising the steps of: 

i) calculating the interpolation coefficients a^, a^,, and a^ respectively corresponding to the original 
image signal components Y k . 1( Y k , Y k+1< and Y k+2 such that, in cases where the original image signal compo- 
nents representing the two picture elements X^ and X,< located before the interpolation point Xp, which is 
located between the picture elements Xk and X k+1 in the original image, are represented by Y M and Y k , the 
original image signal components representing the two picture elements X k+1 and Xk +2 located after the inter- 
polation point Xp are represented by Y k+1 and Y k+2t the third-order spline interpolating function between the 
picture elements Xk and Xk +1 is represented by f k , the third-order spline interpolating function between the pic- 
ture elements X^ and Xk is represented by f k-1 , and the third-order spline interpolating function between the 
picture elements X k+1 and Xk +2 is represented by f k+1 : 

a) the spline interpolating function f k at the picture elements Xk and Xk +1 may satisfy the original image sig- 
nal components Y k and Y k+1 , as represented by Formulas (7) and (8): 

M* k )=Y k (7) 
f k( X k+l)= Y k+l (8) 
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b) the first-order differential coefficient of the spline interpolating function f k at the picture element X k may 
coincide with the first-order differential coefficient of the spline interpolating function f k _! at the picture ele- 
ment X k . as represented by Formula (9): 

c) the first-order differential coefficient of the spline interpolating function f k at the picture element X k+1 may 
coincide with the first-order differential coefficient of the spline interpolating function f k+1 at the picture ele- 
ment X k+1 , as represented by Formula (10): 

V(X k+l )=f k+1 XX k+1 ) 00) 

d) the first-order differential coefficient of the spline interpolating function f k at the picture element X k may 
have an inclination of an arbitrary parameter a with respect to the gradient of the original image signal 
components Y k-1 and Y k+1 representing the picture elements X k ^ and X k+1 , which are located before and 
after the picture element X k , as represented by Formula (11). the arbitrary parameter a being selected pre- 
viously, and determining the sharpness of a secondary image represented by the interpolation image sig- 
nal, 

f k (X k )= a(Y k+1 -Y k ., )/(X k+1 -X k _, ) (11) 

and 

e) the first-order differential coefficient of the spline interpolating function f k at the picture element X k+1 
may have an inclination of said parameter a with respect to the gradient of the original image signal com- 
ponents Y k and Y k+2 representing the picture elements X k and X k+2 . which are located before and after the 
picture element X^i. as represented by Formula (12): 

^XX k+1 )=a(Y k+2 -Y k )/(X k+2 -X k ) (12) 

ii) calculating interpolation coefficients b^, b|<, b k+1 , and b k+2 respectively corresponding to the original image 
signal components Y k _! , Y kr Y k+1 , and Y k+2 such that: 

f) the spline interpolating function f k at the picture elements and X k+1 may slightly deviate from the orig- 
inal image signal components Y k and Y k+1 in accordance with an arbitrary parameter p, as represented by 
Formulas (17) and (18), the arbitrary parameter p being selected previously and determining the sharp- 
ness of the secondary image represented by the interpolation image signal, 

f k (X k )=-0.5p Y k . t +(1 +p)Y k -0.5p Y k+1 (1 7) 

f k<* k + i )=-0.5p Y k+ (1 +p)Y k+1 -0.5PY k+2 (18) 

g) the first-order differential coefficient of the spline interpolating function f k at the picture element X k may 
coincide with the first-order differential coefficient of the spline interpolating function f^ at the picture ele- 
ment X k , as represented by Formula (9): 

V<X k )=f k V(X k ) «(9) 

h) the first-order differential coefficient of the spline interpolating function f k at the picture element X k+1 
may coincide with the first-order differential coefficient of the spline interpolating function f k+1 at the picture 
element , as represented by Formula (10): 

f k *(X k+1 )=f k+1 '(X k+1 ) (10) 

i) the first-order differential coefficient of the spline interpolating function f k at the picture element Xj< may 
coincide with the gradient of the original image signal components Y k ^ and Y k+1 representing the picture 
elements X^ and X k+1 , which are located before and after the picture element X k , as represented by For- 
mula (19): 

fk'(X k )=(Y k+r Y k . 1 )/(X k+1 -X k . 1 ) (19) 
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and 

j) the first-order differential coefficient of the spline interpolating function f k at the picture element X k+1 may 
coincide with the gradient of the original image signal components Y k and Y k+2 representing the picture 
elements X|< and X k+2 , which are located before and after the picture element X^, as represented by For- 
mula (20): 

^k'(X K+1 )=(Y k+2 -Y k )/(X k+2 -X k ) (20) 

iii) calculating a mean value of the set of the interpolation coefficients a^., and b k _ 1( a mean value of the set of 
the interpolation coefficients a k and t^, a mean value of the set of the interpolation coefficients a k+1 and b|< +1 , 
and a mean value of the set of the interpolation coefficients a^ and b k+2 , which sets respectively correspond 
to the original image signal components Y k . 1f Y k , Y k+1 , and Y k+2f the calculated mean values being taken as 
the values of new interpolation coefficients a k . 1 , a k , a k+1l and a k+2 , and 

iv) calculating the interpolated image signal component Y p corresponding to the interpolation point X p in 
accordance with the calculated new interpolation coefficients a k . 1( a k , a k+1 , and a k+2 , and the original image 
signal components Y^ , Y k , Y k+1 , and Y k+2 . 

30. A spline interpolating operation method for an image signal in as defined in Claim 29 wherein the calculated new 
interpolation coefficients a k _i, a k , a k+1 , and a k+2 , which represent the mean values, are respectively represented by 



Formulas (25), (26), (27), and (28): 

a k-1 H(cx+2P+1 )/4}t 3 +{(2a+3p+2)/4}t 2 -{(cc+1 )/4}t-p/4 (25) 

a k ={(-a+6p+7)/4}t 3 +{(a-9p-11)/4}t 2 +P/2+1 (26) 

a k+1 ={(a-6p-7)/4}t 3 +{(-2a+9p+10)/4}t 2 +{(a+1)/4}t-p/4 (27) 

a k+2 ={(a+2p+1)/4}t 3 -{(a+3p+1)/4}t 2 (28) 



in which t, where 0^t<1 , represents the position of the interpolation point X p , the position being taken with respect 
to the picture element X k , that serves as a reference position, and in the direction heading toward the picture ele- 
ment X k+1 , the lattice interval of the original image signal components being set to be equal to 1. 

31. A spline interpolating operation method for an image signal as defined in Claim 29 wherein the original image is a 
radiation image. 

32. A spline interpolating operation method for an image signal as defined in Claim 31 wherein the radiation image has 
been stored on a stimulable phosphor sheet. 

33. A spline interpolating operation method for an image signal as defined in Claim 32 wherein the original image sig- 
nal representing the radiation image is obtained by exposing the stimulable phosphor sheet to stimulating rays, 
which cause the stimulable phosphor sheet to emit light in proportion to the amount of energy stored thereon during 
its exposure to the radiation, and photoelectrically detecting the emitted light. 

34. A spline interpolating operation method for an image signal as defined in Claim 29 wherein the original image has 
been recorded on photographic film. 

35. A spline interpolating operation apparatus for an image signal, wherein an interpolating operation using a third- 
order spline interpolating function is carried out on an original image signal, which represents an original image and 
is made up of a series of original image signal components, and an interpolation image signal is obtained with For- 
mula (6), said interpolation image signal being made up of a series of image signal components, which occur at 
intervals different from those of the original image signal components, 

Y p = a k-1 Y k .i +a k Y k +a k+1 Y k+1 +a k+2 Y k+2 (6) 

in which Y k . 1 , Y k , Y k+1 , and Y k+2 represent the original image signal components representing picture elements X^ 
1? \, X k+1 , and X k+2 in the original image, Y p represents the interpolated image signal component corresponding 
to an interpolation point Xp located between the picture elements Xj< and X k+1 , and a k . 1§ a*, a^, and a^ repre- 
sent the interpolation coefficients, 
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the apparatus comprising: 



1) a storage means for storing information representing the interpolation coefficients a k _ 1t a^ a^, and a^, 
which respectively correspond to the original image signal components Y k _ 1( Y k , Y k+1 , and Y k+2 representing 
the picture elements X k . 1f X^ X k+1 , and X k+2 in the original image, and which are set to be respectively repre- 
sented by Formulas (13), (14), (15), and (16): 

a k . t »-(ay2)t(M) 2 (13) 

a k =(2- a/2)t 3 -(3- a/2)t 2 + 1 (14) 

a k+1 =(a/2-2)t 3 +(3-a)t 2 +(a/2)t (15) 

a k+2 =(a/2)t 2 (M) (16) 

in which a represents the parameter determining the sharpness of a secondary image represented by the 
interpolation image signal, and t, where 0^t<1 , represents the position of the interpolation point Xp, the position 
being taken with respect to the picture element X^ that serves as a reference position, and in the direction 
heading toward the picture element X k+1 , the lattice interval of the original image signal components being set 
to be equal to 1 , 

2) an input means for inputting the arbitrary parameter a, which determines the sharpness of the secondary 
image reproduced from the interpolation image signal, 

3) an interpolation coefficient operation means for calculating the interpolation coefficients a k . 1( a k , a k+1( and 
a k+2 in accordance with the parameter a, the calculation being carried out from the interpolation coefficients, 
which are stored in said storage means, and the parameter a inputted from said input means, and 

4) an interpolating operation means for storing Formula (6) as the operation formula, and calculating the value 
of the interpolated image signal component Y p , which corresponds to the interpolation point Xp, with Formula 
(6) in accordance with the interpolation coefficients a k . 1( a k , a k+1 , and a k + 2 , which have been calculated by said 
interpolation coefficient operation means, and the original image signal components Y k . 1f Y k , Y k+1 , and Y k+2 . 

36. A spline interpolating operation apparatus for an image signal as defined in Claim 35 wherein the original image is 
a radiation image. 



37. A spline interpolating operation apparatus for an image signal as defined in Claim 36 wherein the radiation image 
has been stored on a stimulable phosphor sheet 

38. A spline interpolating operation apparatus for an image signal as defined in Claim 37 wherein the original image 
signal representing the radiation image is obtained by exposing the stimulable phosphor sheet to stimulating rays, 
which cause the stimulable phosphor sheet to emit light in proportion to the amount of energy stored thereon during 
its exposure to the radiation, and photoelectrically detecting the emitted light. 

39. A spline interpolating operation apparatus for an image signal as defined in Claim 35 wherein the original image 
has been recorded on photographic film. 

40. A spline interpolating operation apparatus for an image signal, wherein an interpolating operation using a third- 
order spline interpolating function is carried out on an original image signal, which represents an original image and 
is made up of a series of original image signal components, and an interpolation image signal is obtained with For- 
mula (6), said interpolation image signal being made up of a series of image signal components, which occur at 
intervals different from those of the original image signal components, 

Y p = a k-1 ^ k-1 +a k +a k+1 ^ k+1 +a k+2 ^ k+2 

in which Y k . 1t Y k , Y k+1 , and Y k+2 represent the original image signal components representing picture elements 
1. *k. Xk+1. an d Xk+2 in the original image. Y p represents the interpolated image signal component corresponding 
to an interpolation point Xp located between the picture elements X k and X k+1 , and a k . 1p a k , a^, and a^ repre- 
sent the interpolation coefficients, 
the apparatus comprising: 



1) a storage means for storing information representing the interpolation coefficients a^, a^ a k+1 , and a|< + 2. 
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which respectively correspond to the original image signal components Y k _ 1( Y k , Y k+1 , and Y k+2 representing 
the picture elements X^, X k , X k+1( and X k+2 in the original image, and which are set to be respectively repre- 



sented by Formulas (21), (22). (23), and (24): 

a k . 1 =-(P+1/2)t 3 +(3p/2+1)t 2 -(1/2)t-p/2 (21) 

a k =3(p+1/2)t 3 -(9p/2+5/2)t 2 +p+1 (22) 

a k+1 =-3(p+1/2)t 3 +(9p/2+2)t 2 +(1/2)t-p/2 (23) 

a k+2 =(P+1/2)t 3 -(3p/2+1/2)t 2 (24) 



in which p represents the parameter determining the sharpness of a secondary image represented by the inter- 
polation image signal, and t, where 0^t<1, represents the position of the interpolation point Xp, the position 
being taken with respect to the picture element Xk, that serves as a reference position, and in the direction 
heading toward the picture element X k+1 , the lattice interval of the original image signal components being set 
to be equal to 1 , 

2) an input means for inputting the arbitrary parameter p, which determines the sharpness of the secondary 
image reproduced from the interpolation image signal, 

3) an interpolation coefficient operation means for calculating the interpolation coefficients a^.,. a k , a k+1 , and 
a k+2 in accordance with the parameter p, the calculation being carried out from the interpolation coefficients, 
which are stored in said storage means, and the parameter p inputted from said input means, and 

4) an interpolating operation means for storing Formula (6) as the operation formula, and calculating the value 
of the interpolated image signal component Y p , which corresponds to the interpolation point X pt with Formula 
(6) in accordance with the interpolation coefficients a^, a k( a k+1 , and a k+2 , which have been calculated by said 
interpolation coefficient operation means, and the original image signal components Y k . 1( Y k , Y k+1 , and Y k+2 . 

41. A spline interpolating operation apparatus for an image signal as defined in Claim 40 wherein the original image is 
a radiation image. 

42. A spline interpolating operation apparatus for an image signal as defined in Claim 41 wherein the radiation image 
has been stored on a stimulable phosphor sheet. 

43. A spline interpolating operation apparatus for an image signal as defined in Claim 42 wherein the original image 
signal representing the radiation image is obtained by exposing the stimulable phosphor sheet to stimulating rays, 
which cause the stimulable phosphor sheet to emit light in proportion to the amount of energy stored thereon during 
its exposure to the radiation, and photoelectrically detecting the emitted light. 

44. A spline interpolating operation apparatus for an image signal as defined in Claim 40 wherein the original image 
has been recorded on photographic film. 

45. A spline interpolating operation apparatus for an image signal, wherein an interpolating operation using a third- 
order spline interpolating function is carried out on an original image signal, which represents an original image and 
is made up of a series of original image signal components, and an interpolation image signal is obtained with For- 
mula (6), said interpolation image signal being made up of a series of image signal components, which occur at 
intervals different from those of the original image signal components, 

Y p =a k-1 Y k-1 +3 k Y k +a k+1 Y k+1 +a k+2 Y k+2 (6) 

in which Y k . 1t Y k , Y k+1t and Y k+2 represent the original image signal components representing picture elements X k _ 
1. X k , X k+1 , and X k+2 in the original image, Y p represents the interpolated image signal component corresponding 
to an interpolation point Xp located between the picture elements X k and X k+1 , and a k . 1t a k , a k+1 , and a k+2 repre- 
sent the interpolation coefficients, 
the apparatus comprising: 

1) a storage means for storing information representing the interpolation coefficients a^, a k+1 , and ak+2. 
which respectively correspond to the original image signal components Y k . 1( Y k , Y k+1 , and Y k+2 representing 
the picture elements X^, Xk, X|< +1 , and Xk +2 in the original image, and which are set to be respectively repre- 
sented by Formulas (25), (26). (27). and (28): 
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a k . l =-{(«+2p+1)/4}t 3 +{(2at3p+2)/^t 2 -{(a+1)/4}t-p/4 

a k ={(-a+6p+7)/4}t 3 +{(a-9p-11)/4}t 2 +p/2+1 
a k+l ={(a-6p-7)/4}t 3 +{(-2a+9p+10)/4}t 2 +{(a+1)/4}t-p/4 
a k+2 ={(a+2p+1)/4}t 3 -{(a+3p+1)/4}t 2 



(25) 



(26) 



(27) 



(28) 



in which a and p represent the parameters determining the sharpness of a secondary image represented by 
the interpolation image signal, and t, where 0^t<1 , represents the position of the interpolation point Xp, the 
position being taken with respect to the picture element X^ that serves as a reference position, and in the 
direction heading toward the picture element X k+1 , the lattice interval of the original image signal components 
being set to be equal to 1 , 

2) an input means for inputting the arbitrary parameters a and p, which determine the sharpness of the sec- 
ondary image reproduced from the interpolation image signal, 

3) an interpolation coefficient operation means for calculating the interpolation coefficients a k . 1t a^ a k+1 , and 
a k+2 in accordance with the parameters a and p, the calculation being carried out from the interpolation coef- 
ficients, which are stored in said storage means, and the parameters a and p inputted from said input means, 
and 

4) an interpolating operation means for storing Formula (6) as the operation formula, and calculating the value 
of the interpolated image signal component Y p , which corresponds to the interpolation point Xp, with Formula 
(6) in accordance with the interpolation coefficients a k .-|, a^ a k+1 , and a k+2 , which have been calculated by said 
interpolation coefficient operation means, and the original image signal components Y k . 1# Y k , Y k+1 , and Y k+2 . 

46. A spline interpolating operation apparatus for an image signal as defined in Claim 45 wherein the original image is 
a radiation image. 

47. A spline interpolating operation apparatus for an image signal as defined in Claim 46 wherein the radiation image 
has been stored on a stimulable phosphor sheet. 

48. A spline interpolating operation apparatus for an image signal as defined in Claim 47 wherein the original image 
signal representing the radiation image is obtained by exposing the stimulable phosphor sheet to stimulating rays, 
which cause the stimulable phosphor sheet to emit light in proportion to the amount of energy stored thereon during 
its exposure to the radiation, and photoelectrical^ detecting the emitted light. 

49. A spline interpolating operation apparatus for an image signal as defined in Claim 45 wherein the original image 
has been recorded on photographic film. 
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(54) Interpolating operation method and apparatus for image signals 

(57) An original image signal, which represents an 
original image and is made up of a series of original 
image signal components Yij, is obtained. Correspond- 
ing interpolation coefficients Bij and Cij, which are set 
for each of the original image signal components Yij, in 
two different interpolating functions f and g, represented 
by Formulas (1) and (2), for obtaining two interpolation 
images having different levels of sharpness are linearly 
combined with each other by using Formula (3). A new 
interpolation coefficient Aij is thereby obtained. An inter- 
polating operation is carried out on the original image 
signal components Yij with an interpolating function h, 
which has the new interpolation coefficient Aij and is 
represented by Formula (4). An interpolation image sig- 
nal is thereby obtained, which is made up of a series of 
image signal components occurring at intervals different 
from those of the original image signal components Yij, 
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f=EBij • Yij 
g=LCij • Yij 
Aij=(1-a)Bij+ccCij 
h=EAij • Yij 
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in which i= 1 , 2 and j= 1 , 2 

The coefficient a in Formula (3) is set to be one of real 
numbers including a range smaller than 0 and/or a 
range larger than 1 . 
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